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ABSTRACT 

A stellar wind module has been developed for the PHOENIX stellar atmosphere code for the purpose 
of computing non-LTE, line-blanketed, expanding atmospheric structures and detailed synthetic spectra 
of hot luminous stars with winds. We apply the code to observations of Deneb, for which we report 
the first positive detections of mm and cm emission (obtained using the SCUBA and the VLA), as well 
a strong upper limit on the 850/im flux (using the HHT). The slope of the radio spectrum shows that 
the stellar wind is partially ionized. We report a uniform-disk angular diameter measurement, 0ud = 
2.40 ± 0.06 mas, from the Navy Prototype Optical Interferometer (NPOI). The measured bolometric 
flux and corrected NPOI angular diameter yield an effective temperature of 8600 ±500 K. Least-squares 
comparisons of synthetic spectral energy distributions from 1220 A to 3.6 cm with the observations 
provide estimates for the effective temperature and the mass- loss rate of ~ 8400 ± 100 K and 8±3 x 10~^ 
M0 yr~^, respectively. This range of mass-loss rates is consistent with that derived from high dispersion 
UV spectra when non-LTE metal-line blanketing is considered. We are unable achieve a reasonable fit 
to a typical Ha P-Cygni profile with any model parameters over a reasonable range. This is troubling 
because the Ha profile is the observational basis for Wind Momentum-Luminosity Relationship. 

Subject headings: methods: numerical — radio continuum: stars — stars: atmospheres — 

stars:individual (Deneb) — stars: winds, outflows — techniques: interferometric 

1. INTRODUCTION 

A-type supergiants are the brightest stars at visual wavelengths (up to My — —9) and are therefore among the brightest 
single stars visible in galaxies. For this reason, these stars have been of increasing interest in extragalactic astronomy 
where they show potential as independent distance indicators (Bresolin et al. 2001; Kudritzki et al. 1999, hereafter K99). 
This potential lies in the use of the Wind Momentum-Luminosity Relationship (WMLR) which is derived from Sobolev 
radiation-driven stellar wind theory (K99 and references therein). Testing this relationship and the theory of radiation 
driven winds most critically requires an accurate determination of the stellar mass-loss rates. The WMLR states that 
Mv^ oc where M is the mass-loss rate, Woo is the terminal velocity of the wind, Ri, is a photospheric reference 

radius, and a is a parameter related to the distribution of atomic line strengths for the spectral lines which drive the 
wind. 

The application of the WMLR to nearby galaxies requires a calibration of the relationship using similar stars in the 
Milky Way. To determine mass-loss rates for A-type supergiants, K99 modeled the hydrogen Balmer lines Ha and H7 
for six stars in the Galaxy and M31 using a unified stellar atmosphere model. In their analysis, K99 synthesized the 
observed line profiles by adjusting model values for the mass-loss rate, the velocity law exponent /5, and a line-broadening 
parameter, vt, after first adopting parameters for the effective temperature, Teg, the radius, i?*, the gravity, log (7, and 
the projected stellar rotational velocity. 
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The use of the Ha profile as the sole diagnostic of A-type supergiant wind momenta is problematic because Ha emission 
observed in stars of this type is known to be strongly variable (Rosendhal 1972; Kaufer ct al. 1996). Nevertheless, the 
accuracy of mass-loss rate determinations for A-type supergiants in M31 is claimed to be at the level of 15% (McCarthy 
et al. 1997). However, there has been no consensus on the mass-loss rate for Deneb (a Cyg, HD 197345, spectral type 
A2Ia according to Morgan & Keenan (1973)), the brighest, nearest, and best-studied A-type supergiant, which we might 
expect to serve as an archetype and benchmark. Published estimates of its mass-loss rate range over more than three 
orders of magnitude, from methods which include modeling the line profiles Mg II A2802 (Kunasz & Praderie 1981, 
3.1 X 10-9 - 1.5 X lO"'^ M0 yr-i) and Ha (Kunasz & Morrison 1982; Scuderi et al. 1992, 1.7 ± 0.4 x lO"'' Mg yr'^; 
3.7 ± 0.8 X 10~^ Mq yr~^), analysing low-excitation Fe II lines in the ultraviolet (Hensberge et al. 1982, 1 — 5 x 10"^ 
Mq yr^i), modeling the IR excess (Barlow & Cohen 1977, 6 x 10"'^ M© yr-^), and radio-flux limits (Abbott et al. 1984, 
< 2 X 10"^ Mq yr-i). 

It is not unreasonable to suggest that the degree to which we understand A-type supergiants in general may be tested 
by Deneb. Therefore, it is troubling that the published mass-loss estimates for Deneb vary so widely and this fact casts 
doubt on the what we can infer about the properties of more distant A-type supergiants. The proximity of Deneb allows 
one to more rigorously constrain its physical parameters and apply several difii'erent techniques to estimate its mass-loss 
rate. 

For these reasons, in this paper we investiage Deneb's fundamental properties, and in particular its mass-loss rate. 
To do this we employ a new stellar wind-model atmosphere package developed by Aufdenberg (2000) for use with the 
PHOENIX^ general-purpose stellar and planetary atmosphere code (Hauschildt & Baron 1999, and references therein) and 
bring together new observations of Deneb at visible, infrared, millimeter, and radio wavelengths. 

Our principal motivation for the constructing unified wind models has been that most analyses of early-type supergiants 
still suffer from severe limitations present in the model atmospheres which are employed. Recent systematic analyses 
(Verdugo ct al. 1999) of the atmospheric parameters of A-typc supergiants have employed LTE, line-blankctcd, plane- 
parallel, static model atmospheres. These models do not include the effects (on the temperature structure and the synthetic 
spectrum) of metal line-blanketing and a spherically extended expanding atmosphere. Furthermore, most analyses do 
not allow for departures from LTE, particularly for species heavier than H and He, in computation of the opacity and 
line formation, which are especially important in low-density extended atmospheres. Santolaya-Rey et al. (1997) have 
developed models, employed by K99, which address many of these limitations including the solution of the spheric;al 
expanding transfer equation in the co-moving frame and the non-LTE treatment of hydrogen and helium. Despite 
these important improvements, these models have had limited success in fitting spectral energy distributions (SEDs) 
due to the lack of metal linc-blankcting in both the atmospheric structure and the synthetic spectrum. Fundamental 
stellar parameters derived from line-profile fitting of hydrogen and helium lines alone may not be consistent with the 
SED. Therefore, we believe that fitting both the line spectrum and the SED simultaneously should lead to more robust 
atmospheric parameters for A-type supergiants. 

In §2 we describe the spectrophotometric, spectroscopic, and interferometic observations used here. In §3 limb-darkening 
models and corrections to the angular diamtcr are discussed. The effective temperature, reddening and other parameters 
constrained by the observations are discussed in §4. In §5 we compare the spectrophotometic data from the UV to the 
radio with our synthetic spectral energy distributions. The comparison of our models to portions of Deneb's line spectrum 
is discussed in §6. We summarize our results and conclusions in §7. Details on the model atmosphere construction and 
computation are left for Appendix A. 

2. OBSERVATIONS 

Previous radio observations of Deneb (Drake & Linsky 1989, and references therein) have yielded only upper limits. We 

report positive detections of radiation from Deneb at 1.35 mm with the Submillimetrc Common-User Bolometer Array 
(SCUBA) located on the James Clerk Maxwell Telescope (JCMT) and at 3.6 cm with the Very Large Array (VLA). In 
addition, two of us (KDG & JAS) observed Deneb at 870 /xm with the Heinrich Hertz Telescope (HHT) and we report 
an upper limit on the fltix density at this wavelength. Futhermore, one of us (TEN) measured Deneb's angular diameter, 
previously uncertain at the 30% level, with the Navy Prototype Optical Interferometer (NPOI) and we present an angular 
diameter measurement accurate to 3%. In addition to these new data, we have brought together observations of Deneb 
from Infrared Processing and Analysis Center (IPAC) data archive. The IR, millimeter and radio data are summarized 
in Table 1. 

2.1. SCUBA Observation 

SCUBA was used to observe Deneb on 4 May 1998. Four runs were made in the photometery mode using the 1350 /xm 
filter. Each run consisted of 50 integrations of 18 sec each. Each integration includes both the on-source and off-source 
time, the on-source time is half the value. Deneb was positively detected at 1.35 mm with flux density of 7.8± 1.9 mJy. 

2.2. VLA Observation 

Deneb was observed by 22 antennas of the VLA in the D — >A array configuration on 29 January 1990. The central 
frequency of the observation was 8414 MHz with bandwidth of 50 MHz. The total observation time was 4800 sec, with 
3450 sec on source. Deneb was positively detected at 3.6 cm with a fiux density of 0.23±0.05 mJy. 

^ Also see http://phoenix.physast.uga.edu 
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Table 1 

Infrared. Millimeter. .\nd R.vdio Photometry of Deneb 



Wavelength 


Flux Density 


Observation 


1.25 


609 ± 19 Jy 


COBE/DIRBE band lA " 


2.2 /xm 


268 ± 8 Jy 


COBE/DIRBE band 2 A 
IRTF K band ^ 


2.2 yum 


282 ± 28 Jy 


3.5 fmi 


125 ± 4 Jy 


COBE/DIRBE band 3A ^ 


3.5 /im 


122 ± 12 Jy 


IRTF L band 


4.29 


89.5 ±9.5 Jy 


MSX band Bl " 


4.35 fj,m 


84.5 ±8.0 Jy 


MSX band B2 " 


4.8 


81 ± 8 Jy 


IRTF M band ^ 


4.9 nm 


74.5 ± 2.2Jy 


COBE/DIRBE band 4 


7.76 urn 


33.6 ± 1.7 Jy 


MSX band A " 


10.2 ^im 


22 ± 2 Jy 


IRTF N band ^ 


11.99 urn 


15.6 ± 0.5 Jy 


MSX band B " 


14.55 fim 


11.2 ±0.5 Jy 


MSX band C " 


20.0 


6.7 ±0.7 Jy 


IRTF Q band ^ 


20.68 ^im 


5.38 ± 0.45 Jy 


MSX band D " 


60 /Ltm 


957 ± 287 mJy 


ISOPHOT 


870 ^im 


< 24 mJy 


HHT 


1.35 mm 


7.8 ± 1.9 mJy 


SCUBA 


2 cm 


< 1.0 mJy 


VLA " 


3.6 cm 


0.23 ± 0.05 mJy 


VLA 


6 cm 


< 0.15 mJy 


VLA " 



''http://cobe.gsfc.nasa.gov/cio/ 

''Abbott et al. (1984) 

■^Price et al. (2001); Egan et al. (1999) 

"^This paper 

"Drake & Linsky (1989) 
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2.3. HHT Observation 

Deneb was observed using the 10-m Heinrich Hertz Telescope (e.g., Baars et al. 1999) on Mt. Graham, Arizona on 
22 December 2000. Wc used the MPIfR 19-elemcnt bolometer array which is sensitive at 870 /im = 345 GHz. The 
observations were done during good weather conditions {t225GHz < 0.06). We used source-sky offsets of 100" compared 
to the telescope beam of 22.4". The total on source integration time was 1600 seconds with individual integration times 
of 20 seconds. Mars was observed to provide calibration to physical units. While we did not detect Deneb, we were able 
to put a 3(7 upper limit on this star of 24.0 mJy at 870 /xm. 

2.4. ISO Data Reduction, DIRBE and MSX Data 

Infrared Space Observatory (ISO) observed Deneb at 60 /um with the ISOPHOT instrument (Lemke et al. 1996). This 
data set was downloaded from the ISO archive and reduced using the PIA package (Gabriel et al. 1997). The observations 
consist of three pointings bracketed by observations of the FCSl calibration source. Each pointing was offset by almost 
the full array width from the last pointing and Deneb was observed in the middle pointing. Deneb was centered on the 
middle pixel of the array. The sky levc;l at the position of Deneb was taken as the mean of the sky levels on the middle 
pixel in the two adjacent observations. The flux of Deneb was the difference between the observed flux and the sky level 
so computed. Including the correction for the flux falling outside the central pixel, the 60 /zm flux of Deneb was 957 mJy. 
The uncertainty of this measurement was estimated at 30% (Schulz et al. 1999). 

Infrared photometric data between 1 /im to 5 /im are from the Diffuse Infrared Background Experiment (DIRBE) 
aboard the Cosmic Background Explorer (COBE). These data were obtained from the Astrophysical Data Facility using 
the the DIRBE Point Source Photometry Research Tool. Additional infrared photometric data between 4 /xm and 20 ^,m 
from the Midcourse Space Experiment (MSX) (Price et al. 2001; Egan et al. 1999) were obtained through IPAC. 

2.5. Ultraviolet Data Sources 

Spectrophotometric data of Deneb from the Skylab S-0109 Far-UV Objective Prism Spectrophotometer (Henize et al. 
1979), the Orbiting Astronomical Observatory-2 (OAO-2) spectrometers (Code & Meade 1979), the S2/68 Ultraviolet 
Sky Survey Telescope aboard TDl (Jamar et al. 1976), and the International Ultraviolet Explorer (lUE) are all in close 
agreement regarding the shape and absolute flux level despite being obtained at different epochs (see Table 2). Figure 
1 shows these data over their common interval from 7 /im-i to 4 /im-i (1400 A to 2500 A). The OAO-2 data do show 
systematically higher fluxes relative to the other data sets between 7 /xm~^ and 6 jimT^ . This OAO-2 excess has been 
noted in studies on the UV spectrophotometry of B-type stars (Bohlin & Holm 1984; Aufdenberg et al. 1999) and is 
expected. There is good agreement between lUE and OAO-2 longward of 2500 A until the end of the LWR camera 
beyond 3300 A. The only absolute spectrophotometric data shortward of 1400 A are from lUE. 

Only two large aperture, low dispersion spectra of Deneb were obtained with lUE, one each in the Short Wavelength 
Primary (SWP) and Long Wavelength Redundant (LWR) cameras, exposures SWP 09133 and LWR 07864, respectively. 
We use the lUE data set exclusively in our SED analysis because it has the greatest wavelength coverage, is consistent 
with the other UV data, has probably the best absolute (and relative) flux calibration, and both the long and short 
wavelength spectra were taken on the same date. For the high-dispersion analysis, we have chosen the well-exposed, 
large-aperture spectra SWP 10375 and LWR 07760. The lUE data were obtained from the ESA "lUE Newly Extracted 
Spectra" (INES) system (Rodriguez-Pascual et al. 1999; Gonzalez- Riestra et al. 2000). 

2.6. Additional Spectroscopic and Spectrophotometric Data Sources 



Table 2 

Low-Resolution UV Spectrophotometry of Deneb 



Observatory 


Wavelength Range (A) 


Dates 


# of Scans 


OAO-S^ 


1840-3600 


1969 June 19 


1 


OAO-2 


1400-1800 


1970 May 11 


1 






1970 June 21 


4 


TDl^ 


1360-2740 


1972 March 19 


3 






1972 October 31 




Skylab 


1347-2290 


1973 August 10 


3 


lUE 


1150-2000 


1980 May 26 


1 


lUE 


1850-3300 


1980 May 26 


1 



^ OAO-2 observation dates from Meade (2001). 

''No observation dates available. Dates shown are those of the TDl mis- 
sion. 
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Fig. 1. — Comparison of ultraviolet absolute spectrophotometry of Deneb from Skylab, lUE, OAO-2, and TDl. See text for details. 



6 



Aufdenberg et al. 



Optical spectrophotometric data are from Glushneva et al. (1992) and are in good agreement with data from the 
compilation by Breger (1976) except at 3975 A, where the datum from Glushneva et al. is ~ 50% lower. Two high- 
dispersion {R ~ 20000) optical spectra with wavelength coverage from 3500 A to 5400 A ("blue") and 5775 A to 8630 
A ("red") were obtained on 1 August 1995 with the Heidelberg echelle spectrograph HEROS (Kaufer 1999). Fourier 
Transform Spectroscopy {R ~ 3000) in the K-h&nd is from Wallace & Hinkle (1997). 

2.7. NPOI Observations 

Each interferometric observation consists of a scan made over ten spectral channels spread in increments of equal 
wavcnumber from 850 nm to 650 nm with a mean wavelength of 740 nm. Here we consider only the visibility data from 
the 38-meter East- West baseline, considered the most robust (Nordgren et al. 1999). Deneb was observed 74 times over 
5 nights in 1997 December. Unfortunately, 62 of these observations have no suitable calibration data. Of the remaining 
data, 5 scans from 1997 December 5, 7 scans from 1997 December 6, and 1 scan from 1997 December 16 have suitable 
calibration data taken within one hour of each science scan. These data are listed in Table 3 and plotted in Figure 2. The 
calibration star for each night and its estimated angidar size (Nordgren et al. 1999) are noted in Table 3. For a detailed 
description of the calibration process used for these data see Nordgren et al. (1999). 
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Table 3 
NPOI Deneb Observations 



Observation Time 


Cal. Time" 


Wavelength 


y2 




Spatial Frequency 


(UT) 


(hours) 


(nm) 






(cycles/arcsecond) 




1997 December 5^ 


01:40:05 


0.305 


849.4 


0.5778 ± 


0.0075 


195.5441 


01:40:05 


0.,'!05 


820.9 


0.5572 


0.0074 


202.3308 


01:40:05 


0.305 


793.9 


0.5155 


0.0073 


209.2131 


01:40:05 


0.305 


768.3 


0.5089 


0.0080 


216.1828 


01:40:05 


0.305 


744.2 


0.4803 


0.0072 


223.1807 


01:40:05 


0.305 


722.9 


0.4605 


0.0073 


229.7603 


01:40:05 


0.305 


701.5 


0.4388 


0.0064 


236.7691 


01:40:05 


0.305 


683.1 


0.4165 


0.0068 


243.1457 


01:40:05 


0.305 


664.6 


0.3947 


0.0063 


249.9149 


01:40:05 


0.305 


648.9 


0.3704 


0.0066 


255.9606 


01:52:05 


0.505 


849.4 


0.5914 


0.0078 


192.7878 


01:52:05 


0.505 


820.9 


0.5689 


0.0079 


199.4806 


01:52:05 


0.505 


793.9 


0.5279 


0.0074 


206.2675 


01:52:05 


0.505 


768.3 


0.5210 


0.0083 


213.1379 


01:52:05 


0.505 


744.2 


0.4962 


0.0075 


220.0392 


01:52:05 


0.505 


722.9 


0.4759 


0.0074 


226.5235 


01:52:05 


0.505 


701.5 


0.4544 


0.0068 


233.4355 


01:52:05 


0.505 


683.1 


0.4302 


0.0076 


239.7208 


01:52:05 


0.505 


664.6 


0.4045 


0.0065 


246.3961 


Ul .OZ.UO 


u.ouo 


O'iO . J 


u.oy lu 


U.UUDo 


zOz .004:0 


01:59:49 


0.634 


849.4 


0.5929 


0.0079 


190.9612 


01:59:49 


0.634 


820.9 


0.5655 


0.0081 


197.5905 


01:59:49 


0.634 


793.9 


0.5261 


0.0076 


204.3120 


01:59:49 


0.634 


768.3 


0.5197 


0.0084 


211.1178 


01:59:49 


0.634 


744.2 


0.4928 


0.0078 


217.9575 


01:59:49 


0.634 


722.9 


0.4733 


0.0078 


224.3784 


01:59:49 


0.634 


701.5 


0.4519 


0.0071 


231.2249 


01:59:49 


0.634 


683.1 


0.4297 


0.0078 


237.4519 


01:59:49 


0.634 


664.6 


0.4079 


0.0068 


244.0599 


01:59:49 


0.634 


648.9 


0.3907 


0.0074 


249.9678 


02:12:32 


0.846 


849.4 


0.5999 


0.0086 


187.9221 


02:12:32 


0.846 


820.9 


0.5776 


0.0082 


194.4449 


02:12:32 


0.846 


793.9 


0.5388 


0.0080 


201.0607 


02:12:32 


0.846 


768.3 


0.5376 


0.0089 


207.7591 


02:12:32 


0.846 


744.2 


0.5136 


0.0084 


214.4886 


02:12:32 


0.846 


722.9 


0.4924 


0.0080 


220.8083 


02:12:32 


0.846 


701.5 


0.4722 


0.0074 


227.5415 


02:12:32 


0.846 


683.1 


0.4481 


0.0085 


233.6719 


02:12:32 


0.846 


664.6 


0.4306 


0.0071 


240.1741 


02:12:32 


0.846 


648.9 


0.4116 


0.0076 


245.9880 


1997 December 6" 


01:16:52 


0.087 


849.4 


0.4991 ± 


0.0100 


199.7180 


01:16:52 


0.087 


820.9 


0.4534 


0.0088 


206.6525 


01:16:52 


0.087 


793.9 


0.4377 


0.0090 


213.6804 


01:16:52 


0.087 


768.3 


0.4091 


0.0107 


220.8014 


01:16:52 


0.087 


744.2 


0.3942 


0.0095 


227.9508 


01:16:52 


0.087 


722.9 


0.3724 


0.0091 


234.6640 


01:16:52 


0.087 


701.5 


0.3710 


0.0088 


241.8242 


01:16:52 


0.087 


683.1 


0.3338 


0.0095 


248.3401 


01:16:52 


0.087 


664.6 


0.3138 


0.0088 


255.2507 


01:16:52 


0.087 


648.9 


0.2957 


0.0084 


261.4261 


01:44:35 


0.074 


849.4 


0.6075 


0.0103 


193.6147 


01:44:35 


0.074 


820.9 


0.5023 


0.0075 


200.3382 


01:44:35 


0.074 


793.9 


0.4984 


0.0077 


207.1531 


01:44:35 


0.074 


768.3 


0.4752 


0.0085 


214.0514 


01:44:35 


0.074 


744.2 


0.4401 


0.0081 


220.9846 


01:44:35 


0.074 


722.9 


0.4189 


0.0078 


227.4956 


01:44:35 


0.074 


701.5 


0.4107 


0.0077 


234.4356 


01:44:35 


0.074 


683.1 


0.3948 


0.0074 


240.7516 


01:44:35 


0.074 


664.6 


0.3711 


0.0076 


247.4509 


01:44:35 


0.074 


648.9 


0.3683 


0.0081 


253.4416 


02:14:50 


0.048 


849.4 


0.5821 


0.0094 


187.9029 


02:14:50 


0.048 


820.9 


0.5490 


0.0084 


194.4257 


02:14:50 


0.048 


793.9 


0.5181 


0.0086 


201.0415 


02:14:50 


0.048 


768.3 


0.5153 


0.0093 


207.7393 


02:14:50 


0.048 


744.2 


0.4667 


0.0085 


214.4687 


02:14:50 


0.048 


722.9 


0.4568 


0.0085 


220.7885 


02:14:50 


0.048 


701.5 


0.4567 


0.0092 


227.5210 


02:14:50 


0.048 


683.1 


0.4371 


0.0083 


233.6514 


02:14:50 


0.048 


664.6 


0.4057 


0.0083 


240.1537 


02:14:50 


0.048 


648.9 


0.4228 


0.0095 


245.9638 


02:27:25 


0.054 


849.4 


0.6579 


0.0140 


183.3663 


02:27:25 


0.054 


820.9 


0.6186 


0.0130 


189.7327 


02:27:25 


0.054 


793.9 


0.5967 


0.0127 


196.1848 


02:27:25 


0.054 


768.3 


0.6022 


0.0138 


202.7225 
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A uniform-disk (UD) angular diameter ^ud was derived for each calibrated scan by fitting an Airy pattern 

2 

(1) 

to the visibility data, where kw is the spatial frequency, k — 2n/Xis the wavenumber and w is the projected baseline, 
Old = ^ud/2 is the angular radius, and Ji is the Bessel function of first order. Table 4 lists the least squares best fit 
value of 0UD for each scan and the mean value of 0ud for each night. The model visibility curves corresponding to the 
mean and ±1<t values of ^ud for each night are shown in Figure 2. 

The scans of December 5th are very consistent with each other and with the single scan from December 16th which, 
while noisier, yields a consistent angular diameter. The scans from December 6th, however, show considerable scatter and 
each scan is considerably noiser than the December 5th data, yielding higher values for each of the fits. A histogram 
of all the derived ^ud values is shown in Figure 3. Experience has shown that with well-observed stars the distribution 
of angular diameters follows a Gaussian distribution in the limit of many (~100) scans. The weighted mean angular 
diameter for the 12 scans is 

0^ = 2.40 ± 0.06 mas, (2) 

where the uncertainty represents our estimate for the error in the mean. This value is consistent with the published 
measurements by the CERGA group: ^ud= 2.04 ± 0.45 mas (in the spectral band 600 ± 25 nm) (Koechlin & Rabbia 
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^Time between the science source scan and the cahbration scan used to cahbrate it. 
^CaUbration: 7 Peg, ^UD = 0.421 mas 

"^Calibration: a Lac. Ovd = 0.502 mas 
"^Calibration: ^ Gas, ^ud — 0.264 mas 
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Fig. 2. — NPOI visibility data of Deneb with uniform-disk (UD) fits. For each date, a model visibility curve (solid line) for the mean UD 
diameter, 6uD (see Table 4), is plotted along with the corresponding visibility curves for the ± la error in the mean UD diameter (dotted 
lines), (a) 5 December 1997: 0^ = 2.36 ±0.01 mas (4 scans); (b) 6 December 1997: = 2.41 ±0.17 mas (7 scans); (c) 16 December 1997: 
duo = 2.33 ± 0.02 mas (1 scan). 



Table 4 

Uniform-Disk Angular Diameters 

Observation Time ^ud 
(UT) (mas) 



1997 December 5: 6'ud = 2.364 ± 0.011 
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1997 December 16: 6Iud = 2.325 ± 0.015 
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1985) and 0ud= 2.6 ±0.3 mas (in the spectral range 500-650 nm) (Bonneau et al. 1981). Hanbury Brown's interferometer 
at Narrabri could not observe Deneb due to its high declination. 

3. ANGULAR DIAMETER 

Our model atmospheres (described in detail in Appendix A) predict substantial limb darkening for the Deneb parameters 

determined in the remainder of this paper. This motivates further analysis of the NPOI observations in order to obtain a 
more physically based estimate of the angular diameter. A limb-darkened angular diameter, 6'ld , is normally employed 
in the calculation of the fundamental effective temperature 

1/4 

(3) 

where J-q is the bolometric flux (corrected for interstellar extinction) and a is the Stefan-Boltzmann constant. Unfortu- 
nately, in extended atmospheres the effective temperature is not well defined or is ambiguous because the atmosphere is 
not sufficiently compact to have a unique radius (Baschek et al. 1991). Hence, a value for T^a cannot be separated from 
an adopted value for the radius and the corresponding angular size. 

For well resolved sources (e.g., Betelgeuse, Arcturus), observations can directly yield limb-darkening information from 
the shape of the visibility curve at spatial frequencies beyond the first null. Deneb, however, has yet to be observed at 
suHiciently long baselines (> 100 m) for limb- darkening information to be obtained in this way. Therefore, the details of 
Deneb's center-to-limb variation (CLV) must for now be supplied by theory. 

Limb-darkening coefficients provided in the literature (e.g.. Claret (2000)) are derived from plane-parallel models and are 
not appropriate for Deneb. For example. Figure 4a shows that the model CLV from an expanding atmosphere or "wind" 
model differ significantly from that of a hydrostatic model. Unlike the wind model CLV which falls very gradually from 
40% of the central disk intensity to zero, the hydrostatic model CLV shows a fairly well defined edge. The unambiguous 
signature of the predicted wind model CLV lies in the height of the visibility function beyond the first null (see Figure 
4b). Future measurements should be able to test this prediction. 
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Fig. 3. — Histogram of the uniform-disk angular diameter values derived from each of the 12 scans of Deneb by the NPOI (see Table 4). 
The data are binned to 0.05 mas. The different cross-hatched regions specify the date of each scan. 
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The gradual dedinc of the wind model CLV presents a problem for deriving 6'ld- As noted by Baschek et al. (1991), a 
flat intensity decrease toward the outer disk leaves an "intensity radius" , defined by an arbitrary cutoff of the CLV, e.g. 
at 1%, poorly defined. Such a profile prohibits a simple relation between an "optical depth radius", which corresponds 
to reference optical depth along a ray normal to the atmosphere, and an "intensity radius", which corresponds to the 
impact parameter of a ray at the angle of the intensity cutoff. We choose to define the limb-darkened edge to be the angle 
beyond which the model intensity profile would have an immeasurable effect on the observed visibilities. This is defined 
quantitatively by the following procedure. 

The intensity of the stellar disk as a function of the angular radius, I (a), is extracted from the model atmosphere. 
The coordinate system for the solution of the spherical model atmosphere problem is specified by a set of parallel rays 
tangent to nested spherical sliells (e.g., see Figure 1 in Hummc^r & Rybicki 1971; Mihalas et al. 1975). The rays are 
parameterized by the perpendicular distance p from the line of symmetry though the center of the star. A solution to the 
model atmosphere problem provides the specific intensity 7a at wavelength A at depths z along each ray. The observed 
intensity for a given ray is taken to be the intensity at the intersection of that ray with the outermost shell. The relation 
between the directional cosine of the emergent radiation, /i = cos 6*, and the angular displacement on the sky, a, is 

a=Ci-/i2)i/2. (4) 

The spectral response of the NPOI is approximated by the mean the R- and /-bands (Nordgren et al. 1999). Hence, 
after calculating Ix{a) on a fine wavelength grid, a passband- weighted CLV for each filter, 

Ida) = 7 

is computed, where Si (A) is the response function of the photometric system for the filter i with wavelength limits and 
bi. Response functions for the V-, R- and /-bands are taken from Bessell (1990). 

Following Quirrenbach et al. (1996), the model visibilities, Vi{kw), for each passband are computed by numerically 
integrating 



ii^) = ,5. (5) 



/■aiimb 

27r / Jo{kwa)Ii{a)ada 
Jo 



(6) 



where aumb = Q^ld = is the limb-darkened angular radius and Jo is the Bessel function of zeroth order. The 

squared visibility will be insensitive to the low-intensity tail of the CLV beyond some angle. Therefore, we define aumh to 
be the angle beyond which the CLV affects the visibility at 220 arcsec"^ by less than 1%, approximately the precision of 
the NPOI observations. This corresponds to an intensity cutoff at ^5% of the central disk intensity (see Figure 5a). 

Figure 5b shows that well inside the first null, the shapes of UD and LD visibility curves are essentially identical and 
both fit the NPOI data, but correspond to quantitatively different diameters, 6'ud and 6'ld, respectively. The ratio 
Old/ Ovd= Q!ld/q;ud is the limb-darkening correction. The value of auD parameterizes the UD visibility curve, equation 
(1). The average limb- darkening correction from the R- & /-bands is 15.0%. Therefore, we adopt 

6Ild = 1.15 -eiuD = 2.76 ±0.06 mas (7) 

as the limb-darkened angular diameter of Deneb. 

An accurate and precise limb-darkened angular diameter measurement may yield an effective temperature, via equation 
(3), which is not consistent with the visual color temperature or previously determined values for Teg- from model 
atmosphere fitting. A proper interpretation of this Toff can be made only after an appropriate CLV has been modeled (or 
measured) and the relationships between the intensity and optical depth radii have been established. Three different radii, 
and their corresponding angular diameters, are useful for the comparison of our model atmospheres with the observational 
data: (1) /?*, the input radius to the model atmosphere, L = AirR^aT^, which also corresponds to the radius at the base 
of the wind (see §A.2.2); (2) /?(tross = 2/3), the optical depth radius corresponding a Rosseland optical depth of 2/3; (3) 
/iiimb, the wavelength-dependent radius corresponding to the limb-darkened edge of the atmosphere, the impact parameter 
at the angle anmb- 

An example of the relationship between these radii is provided by one of the best fitting models to Deneb's spectral 
energy distribution (see Table 6, §5.2). For this model the ratio Rumb / R{tross = 2/3) = 1.19. This says that the 
intensity radius derived from the CLV is 19% larger than the optical depth radius in the F-band. Furthermore, this 
model is parameterized with /?*, where /i!<,//?,(rRoss = 2/3) = 0.89, and = 9000 K. Therefore, the effective temperature 

corresponding to Rnmh and the same luminosity is 

Tijr'^ = VR*/R\imb = 7780 K. (8) 

This model effective temperature can be directly compared to T^s from equation 3. However, T^^g* is not the most 
physically relevant temperature and it can be easily shown that a plane-parallel, hydrostatic model with this effective 
temperature is a very poor match to the spectral energy distribution. 

The most useful effective temperature is one which corresponds to the radius at a Rosseland optical depth of tross = 2/3. 
In this model, the V-, R- and /-band optical depth radii all differ by less than 1% from the Rosseland depth. The 
corresponding effective temperature is 

= nVR*/R{TRoss = 2/3) ~ 8490 K. (9) 
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Fig. 4. — (a) Comparison of the /-band ccntcr-to-limb variations (CLVs) from two model atmospfieres and a uniform-disk: uniform-disk 
(dotted line), spliciical, hydrostatic model (dashed line), and spherical hydrostatic plus expanding wind model (solid line). The three CLVs 
have been scaled so that their visibility curves match the NPOI data from 5 December 1997 at 220 arcsec"^. (b) Different CLVs can be most 
easily discriminated beyond the first null; the uniform-disk squared visibilities are predicted to be 60% larger than the wind model squared 
visibilities at a spatial frequency of 700 arcsec"^ in the /-band. The spherical, expanding atmosphere model is parameterized by: T* = 9000 
K at /J* = 1.2 X lO"''^ cm and a surface gravity, logg(/J*) = 1.3, equivalent to T^g = 8543 K at /?(tross = 2/3) and \og g(R) = 1.21, see 
equation 9; a mass-loss rate, M = 5.0 X 10~^ Mq yr"-*^; a terminal wind velocity, Vao 
(see §A.2). The spherical, hydrostatic model is parameterized by: T* = 8500 K at R* 
Teff = 8746 K at /?(tross = 2/3) and logg{R) = 1.35. 



225 kms ; and a velocity law exponent, /3 = 3.0 
= 1.2 X 10^^ cm and logg(/J*) = 1.3, equivalent to 
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Fig. 5. — (a) Model center-to-limb variation in the /-band for an spherical, expanding atmosphere model with parameters: T* = 9000 
K at il* = 1.2 X 10^3 cm and logg(il,) = 1.3, equivalent to Tgff = 8490 K at i?(rRoss = 2/3) and logg(R) = 1.19, see equation 9; M 
= 1.0 X 10~^ Mq yr~^; Vaa = 225 kms~^; and /3 = 3.0. The vertical dashed line, labeled cvld, marks the angular radius beyond which 
Ay ^(220 arcsec""*^) < 1% (see text) . The dotted line, labeled ct\jD, marks the angular radius of a uniform-disk with a visibility curve 
that matches the limb-darkened visibility curve at = 0.5. For this model the limb-darkened radius is 15% larger than the uniform-disk 
radius, (b) Corresponding visibility curves for the model I-band CLVs (dashed, limb-darkened disk; dotted, uniform-disk) over plotted on the 
visibilities measured by NPOI on 1997 December 5. 
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Interestingly, the ratio -Riimb/-R(''Ross = 2/3) is larger than the limb-darkening correction, ai^u/avu- This means that the 
angular diameter corresponding to T^^^ is smaller than the uniform-disk diameter. 

^Ross = ^LD ^^^"^"T ^ ^^^^ = 2.32 ± 0.06 mas (10) 

-Klimb 

With more compact stars, the limb-darkening correction is used to scale up the angular diameter from the uniform-disk 

measurement (^r.oss > ^xjd) in order to bring the fundamental effective temperature in line with the color temperature. 
However, in the case of Dcncb (and likely for other BA-supcrgiants) the angular diameter must be scaled down from 
uniform-disk measurement (^^uoss < ^xjd) for this purpose. In effect, the observed uniform-disk diameter is biased toward 
a larger diameters than it otherwise would be in the absence of the gradual drop-off in the center-to-limb profile predicted 
by our wind models. Longer baseline interferometry of Deneb is needed to test this prediction. 

4. OBSERVATIONAL CONSTRAINTS 

4.1. Effective Temperature and Reddening 

The most fundamental atmospheric parameter of a stellar atmosphere is its effective temperature. As shown above, 
Teff must be carefully defined in the case of an extended atmosphere. Published values for Deneb's TeS range from 7635 
K to 10080 K (Verdugo et al. 1999). Efforts to establish an accurate and precise value of T^s have been ongoing for the 
last thirty years. Constraints on Deneb's fundamental effective temperature have been limited by the large uncertainties 
in its angular diameter and in the foreground extinction which enters into the calculation through the bolometric flux. As 
a result, the value of T^g determined from the bolometric flux and the angular diameter, has been uncertain at a level of 
~ 30%. We can now use the more precise (~ 3%) angular diameter obtained by NPOI to better constrain both T^g and 
the foreground extinction. 

The extinction curve toward Deneb is not well established due to the lack of nearly unreddened A-type supergiants with 
which to perform a pair method extinction curve analysis. As a best guess, we have used an extinction curve derived for 
the well-studied line-of-sight toward the 09 V((f)) star HD 199579 (Fitzpatrick & Massa 1990). Deneb and HD 199579 
lie in the same extinction field, demarcated (84/0.0) by Neckel et al. (1980). The correlation (Jenniskens & Greenberg 
1993) between the linear slope component (c2 = 0.898 for HD 199579, Fitzpatrick & Massa 1990) of the UV extinction 
curve and Ry, the ratio of total to selective extinction, is 

Rv = (4.61 ± 0.22) - (1.89 ± 0.34)c2 (11) 

which yields, 

i?y=2.9±0.5. (12) 
The color excess, E{B — V), toward Deneb is not well constrained in the literature. Published values range from 0.04 
(Weaver & Torres-Dodgen 1995) to 0.10 (Takeda 1994) to 0.15 (Wclty et al. 1994; Abbott et al. 1984). The NPOI angular 

diameter, corrected to the y-band (equation 10), allows us to estimate E{B — V) through the empirical surface brightness 
to intrinsic {V — R) color relation of Barnes et al. (1976). This is done by finding a simultaneous solution to the y-band 
surface brightness 

Fv = 4.2207- 0.1 Vo - 0.51ogeRoss (13) 
and the Barnes-Evans relation in the color range appropriate for Deneb; 

iV = 3.977- 1.390(F-i?)o, -0.17 < (F - i?)o < 0.00 (14) 
These two equations can be expressed in terms of the observed V magnitude, the observed {V — R) color, and E{B — V): 

Fv = 4.2207 -0.l\v-Rv E{B - V) \ - 0.51og6iRoss (15) 

Fv = 3.977 - 1.390 \{V - R) - 0.82E{B - V)] (16) 

using the relationship E{V — R) = Q.82E{B — V) from Barnes et al. (1978). With the observed V and R magnitudes 
(Johnson et al. 1966) and Ry (equation 12) the solution equations (15) and (16) yields 

E{B-V)=Omt°o°ot (17) 
The imccrtainty in E{B — V) is primarily due to the scatter (0.025 mag) in the Barnes-Evans relation, but also to the 
uncertainties in 0ross and Ry- This method of estimating the reddening yields a best estimate for E{B — V); however 
the range of uncertainty still corresponds to the range of values in the literature. 

With best estimates for the angular diameter and reddening we proceed calculate Deneb's effective temperature. The 
UV extinction curve of HD 199579 is used to deredden Deneb's UV spectrum, while the optical data are dereddened using 
the extinction relations from Cardelli et al. (1989). The bolometric flux is calculated by integrating the dereddened UV 
and optical spectrophotometric data from 1150 A to 10725 A. Fluxes at longer and shorter wavelengths provide a negligible 
contribution to the bolometric flux. Spectra SWP08133 and LWP07864 were used for the range 1150 A < A < 3260 A, 
and the optical and near-infrared data from Glushncva et al. (1992) were used for the range 3260 A < A < 10750 A. The 
uncertainty in the bolometric flux is primarily a combination of uncertainties in the extinction, the monochromatic flux 
measurements, and the absolute flux calibrations. Absolute calibration uncertainties are taken to be 10% in the UV and 
5% in the optical, while uncertainties in the angular diameter and extinction are taken from equations (10), (12), and 
(17) respectively. This yields a fundamental effective temperature for Deneb of 

T^R°«« = 8600 ± 500 K. (18) 
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4.2. Photospheric Radius and Gravity 

Taking the Rosseland mean angular diameter of Deneb from equation (10) and the mean Hipparcos parallax, tt = 
1.01 ± 0.57 mas (Ferryman ct al. 1997), yields a radius of i?(rRoss = 2/3) ~ 1.7 x 10^''' cm = 240 Rq. The large 
uncertainty in the parallax causes the distance probability distribution function to have a long tail extending to large 
distances. If the most probable value from the distance probability distribution function is adopted (685 pc, as opposed 
to the mean, 990 pc), the radius is R{thoss = 2/3) ~ 1.2 x 10^^ cm = 172 Rq. 

Humphreys (1978) assigns Deneb to the Cyg 0B7 association, with an adopted true distance modulus of 9.5 (794 
pc), which yields i?(TRoss = 2/3) ~ 1.38 x lo'^^ = 196 Rq. We adopt i?(rRoss = 2/3) = 180 Rq. This radius and 
T^^^ = 8600 K yield a luminosity L = 1.6 x 10^ Lq. For this luminosity, evolutionary tracks of Schaller et al. (1992) and 
(Heger & Langer 2000, considering the effects of rotation and the composition gradient) suggest Deneb's mass lies in the 
range 20 - 25 Mq. Therefore, surface gravities should lie in the range log5'(i?(TR.oss = 2/3)) = 1.2 — 1.4. A full list of the 
stellar parameters discussed in this paper are compiled in Table 5. 

5. deneb's spectral energy distribution 

5.1. Model Parameters 

The most important input parameters which characterize our model atmospheres (for details see Appendix A) are; a 
model effective temperature, T+, at a reference radius, i?*, the gravity, g, at il*^, a terminal velocity, Voo, and beta-law 
exponent, f3, a mass-loss rate, M, a microturbulent velocity, ^t, assumed to be depth independent, and the chemical 
composition. Each of these parameters may be determined to varying degrees of precision by the data. In this section we 
therefore c:on(iiict scnisitivity tc;sts in order to examine how well each parameter is constrained by the observations. 

As outlined by Venn (1995a,b), a number of studies find that metal abundances for Galactic, early A-type supergiants 
are near solar (±0.2 dex). Nitrogen abundances, while enhanced in A-type supergiants relative to unevolved B-stars, are 
consistent within the uncertainties with solar abundances, since the sun is metal-rich relative to the nearby B-type stars 
(see, e.g. Edvardsson et al. 1993). A detailed abundance study of Deneb using LTE plane-parallel model atmospheres 
was recently published by Albayrak (2000). While we would very much like to understand how the accuracy of these 
abundance determinations is influenced by the lack of non-LTE, spherical, and expanding atmosphere effects, this is 
beyond the present scope of this work. In this paper we use solar abundances (Grevesse & Noels 1993) for all models. 

Table 5 

Stellar Parameters for Deneb 
Parameter Symbol Value Reference 



Measured values 



Uniform-disk angular diameter 


0VD 


2.40 ± 0.06 mas 


§2.7, eqn. 2 




Bolometric flux (reddened) 


T 


7.1 ± 0.4 X 10-6 erg cm-^ s"! 


§4.1 




Parallax 


TT 


1.01 ±0.57 mas 


Perryman et al. (1997) 


Derived values 


Limb-darkened angular diameter 


fLD 


2.76 ± 0.06 mas 


§3, eqn. 7 




Rosseland mean angular diameter 










interferometric model 


Sross 


2.32 ± 0.06 mas 


§3, eqn. 10 




SED fit 


^Ross 


2.35 ± 0.01 mas 


§5.2.1 




Color Excess 










Barnes-Evans relation 


E{B-V) 




§4.1, eqn. 17 




SED fit 


E{B - V) 


0.06 ± 0.01 


§5.2.1 




Total/sclcctivc extinction 


Rv 


2.9 ±0.5 


§4.1, eqn. 12 




Bolometric flux (unrcddcncd) 


ro 


9.8 ± 2.3 x 10-^ erg cm-^ s"! 


§4.1 




Effective Temperature 










interferometry 


^Ross 

"2^*1^3 S 


8600 ± 500 K 


§4.1, eqn. 18 




SED x^ fit 


8420 ± lOOK 


§5.2.1, eqn. 22, 


Fig. 6 


Radius 


R(rRof =2/3) 


~ 180 Rq 


§4.2 




Luminosity 


L 


~ 1.6 ±0.4 X 10^ Lq 


§4.2 




Mass 


M 


~ 20 - 25 Mq 


§4.2 




Surface gravity 








R(-rRoss = 2/3) and M 


logs(R(rRo,, =2/3)) 


~ 1.3 


§4.2 




stability, Balmer edge 


logg(il(TRo,, =2/3)) 


1.1 - 1.6 


§5.2.2, Table 7 




Mass-loss Rate 










SED x^ fit 


M 


8±3 X 10-'^ Mq yr-i 


§5.2.1 




spectral lines 


M 


lO"'^ - 10"*^ Mq yr-1 


§6 




Terminal Velocity 


Voo 


~ 225 kms-l 


§6.1, Fig. 19 




VclocitA' Law Exponent 


9 


~ 3 


§5.2.2. eqn. A2, 


, Table 7 


Mit:r()l ui l)ui(Uil \"(l()cil\' 


s/ 


~ 15 kms ^ 


if, ■> -> Table 7 
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5.2. Least-squares Analysis 

To compare quantitatively a set of model spectra to the observed spectrum of Deneb from the UV to the radio, both 
the model and observed spectra have been binned to 20 A in the UV and the model spectra have been binned to 50 
A in the optical and near-IR to match the resolution of the spectrophotometry in this region. At longer wavelengths, 
monochromatic fluxes from model spectra are interpolated at wavelengths corresponding to DIRBE, MSX, ISO, SCUBA, 
and VLA bands. 

The model and the observed SED are compared in each wavelength bin: 



265 / ~" - , , _..\ 2 



2 



El ./e(Ai) - -Fmodel(Ai)|(/^, 1 



where, 



/S(A.) = /e(A0-10"'^^^^(^-^^ (20) 

is the unreddened flux measured at the earth in wavelength bin i, afo{Xi) is the la uncertainty in this value. In bin i, 

the interstellar extinction is A\./Av and the model flux is Fmodei(Ai). The angular diameter, 0* in radians, scales the 
model surface flux (at R^,) to the measured flux. 

0* = ^Ross (21) 

5.2.1. Constraining and M 

A least squares sum was computed for set of 45 models with a range of effective temperatures and mass-loss rates (the 
T^, M grid is shown in Figure 6b). A combination of parameters 9^ (range: 1.5 3.5 mas) and E{B — V) (range: 0.05-0.12 
mag) was found which minimized the for each model. The models, sorted by x^i are listed in Table 6 along with best 
fitting values for 9i, and E{B — V) at each and M. 

Although the best fit model is a good match to the spectrum (see Figures 7 and 8), the reduced chi-square value, 
~ 4.4, is large and the test yields a very low probability of agreement between the observed and theoretical 
distributions indicative of a bad fit. In cases like this Bevington (1969) suggests the use of the F test because: "the 
statistic x^ measures not only the discrepancy between the estimated function [model spectriim] and the parent function 
[actual spectrum] , but also the deviations between the data [observed spectrum] and the parent function [actual spectrum] 
simultaneously." The F test concentrates on the discrepancy between the model spectrum and actual spectrum. The 
probability that model n with F = Xn/Xmin as good a flt to the data as the best fit model is listed in Table 6. From these 
probabilities the la and 2a contours are plotted in Figure 6a. The la region encloses the parameter space T^, ~ 8925 ± 100 
K and M ~ 8 ± 3 x 10~^ Mq yr~^. This effective temperature range is equivalent to 

rp^ss ^ g420 ± 100 K (22) 

after scaling from to 7?(rrjo^^ = 2/3). 

This effective tcnipcaature is in good agreement with the interfcrometric Rosseland fundamental effective temperature 
from equation 18. The least-squares determination of T^^^ has a smaller random uncca'tainty because a correction 
for the interstellar extinction is incorporated into the fitting procedure (with the value for Ry fixed), while the full 
uncertainties in both Ry (equation 11) and E{B — V) (equation 17) from the Barnes-Evans relation enter into the 
interfcrometric value of T^^^. We note that the derived reddening for the best fitting models (Table 6, column 6) yield 
values, E{B — V)~ 0.06 ± 0.01, at the low end of the range of equation 17. This result is consistent with the lower mean 
least-squares T^'^'^ relative to the mean interfcrometric T^^^, since the derived bolometric fiux will be lower when a lower 
interstellar extinction is adopted. In addition, the derived angular diameters for the best fitting models (Table 6, column 
5) correspond (via equation 21) to ^ross = 2.35 ± 0.01 mas, which is in good agreement with interfcrometric value of 6ross 
from equation 10. In short, the derived interferometric and spectrophotometric values for T^^^, Oross, and E{B — V) are 
in good agreement. 

Figure 6b shows the contours of the x^ values from Table 6. Models in the region (8500 K, 2x10^** M0 yr^^) poorly 
fit the observed spectrum because the model UV fluxes are too small, particularly right below the Balmer jump and 
below 1600 A, and the model millimeter and radio fluxes are too small by roughly one order of magnitude. Models with 
much higher mass-loss rates and similar effective temperatures (8500 K, 2x10^^ Mq yr~^) provide a closer match to the 
millimeter and radio fluxes, but still come up short, and the model UV fluxes are even smaller. Models in the region (9250 
K, 2xl0~^ Mq yr~^) poorly flt the observed spectrum because the model UV fluxes below 1600 A are too large, and 
the predicted millimeter and radio fluxes are too small by an order of magnitude. For a given T^, larger M values have 
two effects on the SED relative to lower M values: (1) the UV continuum is suppressed, and (2) the radio continuum is 
enhanced. The hottest, highest mass-loss rate models explored (eg., 9250 K, 3xl0~^ Mq jt ) produce millimeter and 
radio fluxes in excess of those observed, and the model UV continuum, while suppressed relative to lower M models, is 
still too large below 1600 A. The best fitting models are slightly cooler, to match the UV, and have slightly lower mass-loss 
rates, to match the millimeter and radio fluxes. 
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Table 6 

SED (UV TO Radio) Fitting Results" 





Model 




Fit Results 




T* (K) 


M (M0 yr ) 


X 


P{F test) 


Oi, (mas) 


E{B — V) 


8875 


l.OOE-06 


1.14E-I-03 


5.00E-01 


2.10 


0.05 


9000 


l.OOE-06 


1.17E+03 


4.13E-01 


2.08 


0.06 


8875 


5.00E-07 


1.21E+03 


3.18E-01 


2.10 


0.05 


9000 


5.00E-07 


1.23E-)-03 


2.57E-01 


2.08 


0.06 


8750 


3.00E-07 


1.24E4-03 


2.44E-01 


2.16 


0.05 


8875 


3.00E-07 


1.30E-I-03 


1.44E-01 


2.16 


0.07 


9125 


l.OOE-06 


1.31E+03 


1.26E-01 


2.08 


0.08 


9000 


2.00E-06 


1.32E+03 


1.21E-01 


2.06 


0.05 


8875 


2.00E-06 


1.34E+03 


9.28E-02 


2.10 


0.05 


9125 


5.00E-07 


1.37E4-03 


7.08E-02 


2.10 


0.09 


8750 


l.OOE-07 


1.38E-I-03 


5.97E-02 


2.20 


0.06 


9000 


3.00E-07 


1.43E+03 


3.22E-02 


2.16 


0.09 


9125 


3.00E-07 


1.45E+03 


2.72E-02 


2.10 


0.09 


9125 


2.00E-06 


1.45E+03 


2.50E-02 


2.06 


0.07 


8750 


l.OOE-06 


1.47E-I-03 


1.99E-02 


2.16 


0.05 


8750 


5.00E-08 


1.47E4-03 


1.98E-02 


2.26 


0.08 


8750 


3.00E-08 


1.49E+03 


1.46E-02 


2.18 


0.05 


9000 


3.00E-06 


1.52E+03 


9.91E-03 


2.04 


0.05 


8500 


l.OOE-07 


1.52E+03 




2 28 


0.05 


9000 


2.00E-08 


1.54E+03 


7.58E-03 


2.18 


0.09 


9250 


l.OOE-06 


1.55E+03 


6.71E-03 


2.08 


0.10 


8500 


5.00E-08 


1.55E+03 


6.38E-03 


2.28 


0.05 


8750 


2.00E-08 


1.57E+03 


4.53E-03 


2.18 


0.05 


8750 


5.00E-07 


1.59E+03 


3.53E-03 


2.16 


0.05 


9250 


5.00E-07 


1.61E+03 


2.75E-03 


2.10 


0.11 


9000 


l.OOE-07 


1.62E+03 


2.23E-03 


2.20 


0.10 


9000 


3.00E-08 


1.63E+03 


1.89E-03 


2.20 


0.10 


9125 


3.00E-06 


1.64E+03 


1.65E-03 


2.02 


0.06 


8750 


2.00E-06 


1.64E+03 


1.65E-03 


2.14 


0.05 


8875 


3.00E-06 


1.65E+03 


1.32E-03 


2.10 


0.05 


9250 


2.00E-06 


1.67E-t-03 


9.48E-04 


2.06 


0.09 


9250 


3.00E-07 


1.70E4-03 


6.61E-04 


2.10 


0.11 


9000 


5.00E-08 


1.72E-F03 


4.52E-04 


2.22 


0.11 


9250 


3.00E-06 


1.86E+03 


3.70E-05 


2.00 


0.07 


8500 


3.00E-07 


1.90E+03 


1.74E-05 


2.28 


0.05 


8500 


5.00E-07 


1.95E+03 


8.16E-06 


2.26 


0.05 


8500 


3.00E-08 


2.00E-I-03 


2.98E-06 


2.30 


0.05 


8750 


3.00E-06 


2.17E4-03 


1.19E-07 


2.12 


0.05 


9250 


l.OOE-07 


2.69E+03 


0.00 


2.14 


0.12 


9250 


2.00E-08 


2.71E+03 


0.00 


2.14 


0.12 


9250 


3.00E-08 


2.81E-t-03 


0.00 


2.14 


0.12 


8500 


3.00E-06 


2.93E-F03 


0.00 


2.22 


0.05 


8500 


2.00E-06 


2.94E-F03 


0.00 


2.22 


0.05 


9250 


5.00E-08 


3.11E+03 


0.00 


2.12 


0.12 


8500 


l.OOE-06 


3.25E+03 


0.00 


2.24 


0.05 


8500 


2.00E-08 


3.36E-t-03 


0.00 


2.28 


0.05 



^Columns (1-2): With the exception of parameters and M. all models listed 
share these parameters: ii* = 1.2x10^^ cm, \ogg{R^) — 1.3, ^ — 15 kms~^, v^o 
= 225 kms^^, /3 = 3.0. Column (3): values are from equation (19) for 265 
wavelength bins from 1220 A to 3.6 cm. Column (4): The probability from the 
F test that the value for each model is statistically equivalent to the best fit 
model. Column (5): The best fit angular diameter value for each T*. , M combination. 
Column (6): The best fit color excess for each T^, Mcombination. i?v=2.9 is fixed. 
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1 I . . . , . . . , . . . , . . . , 

8400 8600 8800 9000 9200 

Model Effective Temperature, T, [K] 



Fig. 6. — (a) Regions (Icr and 2a) in the (T,, M) parameter plane containing the best fit model parameters to the observed spectral energy 
distribution (SED) of Deneb. The probabilities axe listed in Table 6. (b) Contour map of in the (T*,M) plane of the model fits to the 
observed SED. The grid of model locations in the parameter plane is marked by x. 
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Fig. 7. — Top Panel: Best least squares fit (solid line) to the UV to radio spectrum compared to the UV-optical spectrum of Deneb (error 
bars), which is binned to 20 A in the UV and to 50 A in the optical. The model has been scaled with the angular diameter 9* (see text) and 
the data have been corrected for interstellar extinction using the reddening parameters shown. T* = 8875 K is equivalent to T^^^ = 8370 K. 
Bottom Panel: Percentage difference, (model - data)/data, plotted as a function of wavenumber. 
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Fig. 8. — Top Panel: Best least squares fit (solid line) to the UV to radio spectrum compared to IR, mm, radio spectrum of Deneb (error 
bars) from Table 1. The model has been scaled with the angular diameter 9^,. = 8875 K is equivalent to T^'^" = 8370 K. Bottom Panel: 
Percentage difference, (model - data)/data, plotted as a function of wavelength. 
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5.2.2. Constraining Other Parameters 

Our preliminary fits to Deneb's SED revealed that model parameters and M have the largest effects on the model 
SED. As a result, the parameters i?*, log^f, Voo, /?, and were fixed in the analysis discussed above. The most straight- 
forward parameter to constrain is the terminal velocity. We find Vx, ~ 225 kms~^ by fitting the blue edge of the Mg II h 
resonance line profile (see Figure 19). Parameters /?, ^t, and loggr can be constrained with reference to least squares fits to 
the observed SED (see Table 7). Parameters f3 and are constrained due to their significant effect on the ultraviolet SED 
(see Figure 9). Despite the uncertainty in the reddening, /3 values 2.0 and 4.0 can be rejected with confidence. Similarly, 
a low microturbulence, = 5kms~^, can be rejected. 

We find that adjusting the model radius iJ* by ±25% has no significant effect on the model SED. However, as shown 
below in §6.2, the radius can affect synthetic P-Cygni profiles because, via the continuity equation (Al), the parameter 
i?*, for a fixed mass- loss rate, sets the gas density at the base of the wind. 

The model with log g(i?^) = 1.7 can be rejected at the 3(t level due to this parameter's strong effect on the Balmer 
jump. From the requirement of the hydrostatic equilibrium (equations A7, A8, and A9) in the deepest layers (tross — 100) 
of the model atmosphere we can assign a lower limit of to the gravity of logg{R^,) > 1.2. Models with lower values of 
log (7(i?Vf) fail to have positive effective gravities in the deepest layers of the atmosphere, thus making a hydrostatic solution 
impossible. The radiative acceleration is quite significant in the hydrostatic layers. For example, in the log g(_R^) = 1.7 
model the effective gravity, log^eff) drops as low as 1.43 in the hydrostatic layers, while for the logg(i?*) = 1.3 model, the 
minimum log^eff is 0.4. The surface gravity constraint is then 1.2 < loggi(i?*) < 1.7. This range of values can be scaled 
to the reference depth i?(rRoss = 2/3) via 

logg{R{Tnoss = 2/3)) = log g{R^) + 2 log {R^/R{moss = 2/3)) (23) 

yielding 1.1 < log(?(i?(TRoss = 2/3)) < 1.6 for the same ratio of of radn (0.89) as equation 9. 

5.3. Deneb's Continuurn, Beyond 1 ^im 

Figure 10 shows the IR, millimeter, and radio observations from Table 1 compared with the synthetic continua from 
a set of T^, = 9000 K models listed in Table 6. The synthetic continua are scaled with 9^, = 2.08 mas, consistent with 
the NPOI diameter. Models with mass-loss rates M > 3 x 10~^ Mq yr~^ can be ruled out by the 3a upper limit at 
870 /xm from the HHT. Except for the DIRBE band lA, observations at shorter wavelengths are consistent with this. 
The detections at 1.35 mm and 3.6 cm, with SCUBA and the VLA respectively, rule out models with mass-loss rates M 
< 5 X 10"'' M0 yr~^. These results are consistent with the least squares results from the full UV to radio spectrum. 

The synthetic continua shown in Figure 10 were computed from spherical expanding atmosphere models whose temper- 
ature structures are consistent with radiative equilibrium and where non-LTE line-blanketing has been self-consistently 
included. This is important because at wavelengths beyond 1 mm the synthetic continua are quite sensitive to the degree 
of line-blanketing employed in the construction of the model atmosphere. This is illustrated in Figure 11, which shows 
the results of a test where five wind models, which differ in the degree of line-blanketing, but have otherwise identical 
parameters, were computed. As easily seen with reference to the hydrostatic model spectrum, the expanding atmosphere 
spectra predict significantly larger fluxes. The effects of line-blanketing and non-LTE are no less significant. For example, 
at 3 cm, the continuum opacity wind model produces ^10 times more fiux than our most sophisticated non-LTE line- 
blanketed wind model, and the LTE line-blanketed wind model produces ^1/2 as much flux. The LTE continuum-only 
model spectrum shows a turnover at ~ 30 cm due to truncation at the outer boundary of the model, 200 i?*. A similar 
model with an other boundary of 400 i?* does not show a turnover. The line-blanketed wind models, which are not fully 
ionized to the outer boundary, are not sensitive to this truncation. 

There are two important factors at play here. First, the incorporation of metal line-blanketing in the radiative equilib- 
rium solution cools the outermost layers of the atmosphere (classic surface cooling) , especially at depths where the radio 
continuum forms. In LTE, lower electron temperatures mean lower ionization, and therefore less thermal radio emission 



Table 7 

Least Squares Results for Parameters ii*, logg, f3, ^t'^ 





Model 








Fit Results 




logs(iJ.) 


R* 




it 




P{F tost) 


e. 


E{B - V) 


1.3 


1.50e-)-13 


3.0 


15 


1.15e+03 


5.00e-01 


2.080 


0.060 


1.3 


1.20e-M3 


3.0 


15 


1.17e+03 


4.32e-01 


2.080 


0.060 


1.3 


9.00e-M2 


3.0 


15 


1.26e+03 


2.25e-01 


2.080 


0.060 


1.3 


1. 200+13 


3.0 


25 


1.31C+03 


1.36O-01 


2.060 


0.050 


1.5 


1. 200+13 


3.0 


15 


1.46O+03 


2.660-02 


2.060 


0.050 


1.3 


1.20e+13 


2.0 


15 


1.69e+03 


8.17e-04 


2.260 


0.120 


1.7 


1.20e+13 


3.0 


15 


1.87e+03 


3.63e-05 


2.040 


0.050 


1.3 


1.20e+13 


3.0 


5 


2.13e+03 


2.98e-07 


2.120 


0.090 


1.3 


1.20e+13 


4.0 


15 


2.23e+03 


5.96e-08 


2.020 


0.050 



"Columns (1-4): The modols listod here have the fixed values T, = 9000 K and M = 
lO"** M0 yr~^ Columns (5-8) as described in Table 6. 
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Fig. 9.- Synthetic UV continua from three sets (from top to bottom) of three models: surface gravity, logg = 1.3 (solid), 1.5 (dotted), 1.7 
(dashed); the microturbulent velocity, = 5 km s~^ (dotted), 15 kms~^(solid), 25 km (dashed); and the velocity law exponent, (3 = 2.0 
(dotted), 3.0 (solid), 4.0 (dashed). These synthetic spectra have been binned to 20 A resolution to show the gross features of the spectrum. 

All models have the following parameters, unless otherwise noted: = 9000 K , log(y(R*) = 1.3, M = 10~^ M,:) yr^^, ,3 = 3.0, = 15 
kms^^. All models have been scaled by = 2.08 mas. Constants of 1.25 X 10"'^ and 2.5 X lO"'-* erg s ^ cm ^ A ^ have been added to the 
top two sets of models for clarity. The effects of log g can be seen on both sides of the Balmcr discontinuity. The effects of varying from 5 
to 25 kms~^ are apparent below 3000 A (> 3.5 /im~^) where the flux drops by roughly a factor of two. The effects of varying /3 from 2.0 to 
4.0 are significant below the Balmer jump where the flux drops by roughly 50%. 
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Fig. 10. — Infrared, millimeter and radio photometry of Deneb from Table 1. Synthetic continua for models listed in Table 6 with = 
9000 K and 6 different mass-loss rates are shown for comparison. The inset shows a close up of the COBE/DIRBE and MSX photometry. 
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Fig. 11. - Comparison of 4 wind models with identical parameters (T, = 9000 K, T^""" = 8600 K, M = 10"*' Mq yr"\ i?* = 1.2el3 
cm, \ogg{Ri,) = 1.3, ^ = 15 kms"-'-, Voo = 225 kms"^, ,i3 = 3.0). sliowing the importance of line-blanketing on the predicted radio flux. A 
hydrostatic model (T* = 8600 K, ij* = 1.2el3 cm, logg(iJ*) = 1.3, ^ = 15 kms~^) is shown for reference. These synthetic spectra have been 
scaled with Oi, = 2.08 mas. The HHT, SCUBA, and VLA data are shown for comparison. 
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is expected in the presence of line-blanketing relative to the continuum-only opacity case. Second, photoionization is en- 
hanced in non-LTE. Non-LTE photoionization affects the ionization state of the outer layers because they are radiatively 
coupled to the hotter radiation field at depth. In non-LTE, the ionization fraction in the outermost layers is enhanced 
relative to the LTE line-blanketed case, which leads to more thermal radio emission. 

The sensitivity of the radio continuum to the degree of ionization in the extended stellar envelope makes it clear that 
the application of a fully ionized, uniform, spherically symmetric mass flow model commonly employed in mass-loss rate 
estimates for O- and B-type supergiants (Scuderi et al. 1998; Drake & Linsky 1989) won't work for A-type supergiants. 
Furthermore, in such models the radio spectrum follows the power-law Sp oc u"^/^, but Deneb exhibits a spectral slope 
between 1.35 mm and 3.6 cm consistent with Si, oc j/i-26±o.oi^ g^j^j therefore inconsistent with a fully ionized wind. 
Simon et al. (1983) studied thermal radio emission from partially ionized circumstellar environments and winds. They 
found that the slope of radio spectrum varies from the fully ionized, optically thin case, {S^ oc v^-^), to the optically 
thick, Rayleigh- Jeans case, (5^ oc v^). Deneb is an intermediate case. The steepness of a radio spectrum depends on the 
degree of ionization and radial extent of the ionized envelope. We find the same result, where the degree of ionization is 
controlled by line cooling and photoionization. The optically thick case is represented by a geometrically thin, hydrostatic 
model (see Figure 11). 

As demonstrated above, for a given efi^ective temperature, the degree of line-blanketing, photoionization, and the mass- 
loss rate have the largest effects on the synthetic spectrum beyond 1 ^m. The synthetic spectrum is also aff'ected to a 
much smaller degree by other model parameters. Figure 12 shows the sensitivity of the synthetic spectrum, relative to a 
good fitting reference model, to changes in /?, T*, i?^, ^j, and \ogg. Changing each of these parameters has systematic 
effects on the synthetic radio continua. Higher (lower) model temperatures increase (decrease) the ionization in the wind 
and produce higher (lower) fluxes. Changing by 250 K has little effect below 1 mm. Larger (smaller) values of R^, 
yield lower (higher) wind densities for a fixed mass-loss rate and as a result yields slightly lower (higher) fluxes. Higher 
log 5 values yield slightly lower fluxes and varying has no significant effect on the fluxes. 
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Fig. 12. — Synthetic IR, millimeter, radio coiitiima from five sets (from top to bottom): the velocity law exponent, (3 = 2.0 (dotted), 
3.0 (solid), 4.0 (dashed); the model tcmpcratm'c, T^, = 8750 K (dotted), 9000 K (solid), 9250 K (dashed); the model radius, _R* = 130 
Rq (dotted), 172 Rq (solid), 215 Rq (dashed); microturbuleiice, = 5 kms^-*-, 15 kms^-*-, 25 kms^-"-; and the surface gravity, logg = 
1.3 (solid), 1.5 (dotted), 1.7 (dashed). All models have the following parameters, unless noted otherwise: T* = 9000 K , log(/(_R*) = 1.3, 
M = 10~^ Mq yr"'^, /3 = 3.0, £,t = 15 kms~^. The models are scaled using (?* = 2.08 mas. Each set of three models is offset by a factor of 
10 for clarity. The DIRBE, MSX, ISO, SCUBA, and VLA data are shown for comparison. 
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The velocity law parameter, f3, has a very interesting effect on the synthetic fluxes. At wavelengths beyond 0.5 mm, 
larger values of /3 result in lower fluxes, while shortward of 0.5 mm the opposite is true. This behavior in the synthetic 
radio spectra is a result of the run of electron pressure with depth in the wind as a function of /3. At shorter wavelengths, 
where the opacity is lower, the continuum forms in the inner wind 3i?j,) where the gas and electron densities are 
regulated primarily by the velocity law through the continuity equation. As a result, models with larger /? values have 
higher densities at this depth and thus higher fluxes. Further out in the wind (~ lOi?*), where the 3 cm continuum 
forms, non-LTE effects are more important and the situation is reversed. At these depths the electron temperature and 
gas density structures are essentially identical for all (3 values, but the electron densities are higher for models with the 
smaller f3 values. This happens because the models with smaller f3 values have warmer temperature structures at depth 
(sec, for example, the Balmer continuum as a function of /3 in Figure 9) and the resulting radiation field from these layers 
ionizes the wind more substantially than models with larger (3 values. For example, the (3 = 2.0 model has the majority 
of hydrogen ionized out to a radius of i? = 1640 Rq, while for the (3 = 4.0 model, the majority of hydrogren is ionized 
out to only R = 610 Rq. 

In summary, of all the model parameters, M has the largest effect on the synthetic continuum beyond 1 /um. Non-LTE 
photoionization, high Tj, values, and small /3 values all increase model fluxes at centimeter wavelengths for a fixed M. 

6. deneb's line spectrum 
6.1. The Ultraviolet Line Spectrum 

Deneb's rich UV line spectrum, with its multitude of iron peak absorption lines, is a good diagnostic for the column 

density in the expanding atmosphere and therefore should be a good diagnostic for the mass-loss rate. Modeling the 
detailed UV spectrum provides an independent mass-loss rate estimate from that of the IR, mm, radio continuum because 
the latter is based on the integrated emissivity of the extended atmosphere, while the former is based on the line-of-sight 
column density toward the stellar photosphere. 

Let us examine a portion of an lUE spectrum relative to a hydrostatic model synthetic spectrum. Figure 13 shows that 
the synthetic spectrum reproduces the overall appearance of the observed spectrum quite well. However, in the observed 
spectrum the stronger lines are considerably less saturated, broader, and blended relative to the model spectrum. This is 
almost certainly due to these lines forming in the wind and sampling its velocity field, as has been noted by others (e.g. 
Praderie et al. 1980). The 1640 A to 1670 A region is chosen because here the synthetic spectra are quite sensitive to 
the value of M and this spectral region is dominated by Fe II absorption lines. In all comparisons to the lUE data, the 
models are scaled with an angular diameter of 6** — 2.08 mas and the data have been corrected for interstellar extinction 
(curve for HD 199579, E{B - V) = 0.06, Ry = 2.9). 

Now we examine the synthetic UV spectra from the same models we compared to the data beyond 1 /zm (sec Figure 
10). Figure 14 shows a comparison of six synthetic spectra with LTE metal- line blanketing and mass- loss rates from 
2 X 10"^ Mq yr-i to 3 x 10"*^ Mq yr'^ to the lUE spectrum SWP 10375. The lowest mass-loss rate model, M = 
2 X 10~^ M0 yr~^, appears to be most consistent with the windows of low opacity in the observed spectrum: 1645 A, 
1656 A, 1660 A, and 1662 A. Spectra from these LTE models with higher mass-loss rates (higher wind column densities) 
show a very distorted spectrum, inconsistent with the observations. In these higher M models, the line spectrum forms 
at larger radii, where the velocity field distorts and washes out the spectrum. However, even the M = 2 x 10~^ Mq yr~^ 
model cannot provide as good a match to the observed spectrum as the hydrostatic model (see Figure 13). These results 
are very much at odds with the IR, mm, radio continuum anaylsis which indicates M > 10~^ Mq yr~^. 

A solution to this inconsistency between the radio and UV mass-loss rate estimates is to use a model which treats 
Fe I-III in non-LTE (see §A.l for details on the construction of the model ions and the atomic data). This is critically 
important. Figure 15 shows that the temperature structures for the LTE and non-LTE Fe i-iii models are quite similiar, 
as expected since the synthetic radio spectra from the two models are in close agreement (sec Figme 11). However, Figure 
16 shows the dramatic non-LTE effect of the radiation field on the ionization structure of Fe in the wind: Fe"*" is ionized 
to Fe"'""'", significantly dropping the Fe"'" column density in the wind relative to the LTE case. Most of the strong Fe II 
transitions between 1640 A and 1670 A have lower levels between 0.2 cV and 1 eV above the ground state, in terms 
a^'F and a"'D. These levels are seriously under-populated relative to LTE as shown in Figure 17. The resultant non-LTE 
synthetic spectrum (Figure 18) is a much improved match to the observed spectrum with a mass-loss rate consistent with 
the radio continuum result. 

The bottom panel of Figure 18 shows the comparison of the non-LTE model spectrum and the lUE data. The overall 
match to the 1640 A to 1670 A region is generally good, but several strong lines do not appear in the model spectrum: 
e.g., 1645.015 A, 1648.403 A, and 1667.913 A. These lines do appear in the hydrostatic model (see Figure 13) and are 
therefore in the LTE line list (Kurucz & Bell 1995), but are missing or very weak in the non-LTE spectrum. Specifically, 
consulting Kurucz (1994) we find that the 1645.015 A line (d^ a^F9/2 (3H)4p z^Iii/2) is missing from the synthetic 
spectrum because the upper level (62662 cm~^) is missing from the model atom. The 1648.403 A line ((5D)4s 
- (3P)4p y'*P3/2) is an intercombination line in LS-coupling, and its predicted gf-valne falls below the \og{gf) > —3.0 
ciitoff for our Fe II model atom. The data suggest the true g/-value is larger, which is not surprising given the large 
uncertaintly likely for a theoretical gf-value from a spin-forbidden transition. The 1667.913 A line (d'' a^Pi/2 - (lD)4p 
w^Pi/2) by comparison has a more highly excited lower level (18882.702 cm~^ versus 1872 cm~^ to 8846 cm~^ for the a*F 
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Hydrostatic model: T. = 8500 K, Log g = 1 .3, = 1 5 km s"' 
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Fig. 13. — Section of the high-dispersion lUE spectrum SWP 10375 (bottom panel) compared with a synthetic spectrum (top panel) from 
a hydrostatic model. There is good correspondance between the strong absorption features in the model and the data, however the synthetic 
line profiles are much too saturated relative to the observed spectrum. The velocity field present in Deneb's atmosphere desaturates and 
blends the strong absorption lines. 
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Fig. 14. - Section of the high-dispersion lUE spectrum SWP 10375 (bottom panel) compared to a set of synthetic spectra (top panel) 
from models with different mass-loss rates. This spectral region is dominated by Fe II lines, however in these models all species except H 1 
and He I-U are treated in LTE. 
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and a^D terms) and may be quenched by the Boltzmann factor, and possibly by inaccurate radiative and coUisional rates. 
In constrast, the strong hnes which are well matched by the synthetic spectrum (e.g., 1643.578 A, 1650.704 A, 1663.222 
A) have upper levels with highly purc-LS eigenvectors (weak intermediate coupling) which arc expected to yield the most 
reliable transition strengths. The complexity of the high-dispersion UV spectrum and the limitations of the atomic data 
do not allow to us to better constrain other model parameters, such as /3, from this comparison. 

A non-LTE treatment also provides a better fit to the Mg II h and k resonance lines and their satellites as shown in 
Figure 19. The non-LTE synthetic spectrum more accurately reproduces the observed spectrum relative to the LTE case 
in two respects: (1) the reduced strength of the h and k P-Cygni emission components and. (2) the lower velocity shifts 
of the Mg II satellite lines at -450 kms^^ and +325 kms^^. As in the case of Fe II, enhanced ionization in the non-LTE 
case reduces the Mg II column density such that the unsaturated satellite lines form deeper in the wind at lower velocities. 
The models in Figure 19 are parameterized with v^o = 225 kms~^ (maximum model velocity, Vmax = 222 kms~^), 
and provide a good match to the blue edge of the h and k P-Cygni absorption troughs. Because of their strength, the 
resonance lines probe the highest velocities in the wind and thus fix the Voo parameter. 

6.2. Visual-IR Line Spectrum 
6.2.1. The Ha line 

At optical wavelengths, the most conspiciious indicator of mass-loss in early-type supcrgiants is the P-Cygni character 
of the Ha line. Like all B A- type supergiants, Deneb's Ha profile is variable (Kaufer et al. 1996). Deneb's Ha profile 
distinguishes itself from most supergiants in that it lacks broad emission wings normally attributed with electron scattering. 
Despite being the protype of the class, Deneb's Ha profile has only been modeled in detail twice to our knowledge. 

Reasonable fits to several of Deneb's Ha profiles have been achieved by Kunasz & Morrison (1982) and Scuderi et al. 
(1992). Unfortunately, the model parameters and assumptions used in these works are inconsistent with the present 
observational data. In the former work, the best fits were achieved by adopting terminal velocities of 120 kms~^ and 140 
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Fic;. 15. — Comparison of the velocity (left ordinate) and electron temperature (right ordinate) structures for two models with and without 
non-LTE metal line-blanketing of species Cl-lll, Nl-Ill, Ol-Ill, Mgll, Fel-Ill. Both models have full LTE line-blanketing for those species not 
treated in non-LTE. The models have the following parameters: = 9000 K , logg(/?*) = 1.3, M = V)-^ Mq yr"\ /? = 3.0, it = 15 
kms-i. The non-LTE model has a slightly warmer temperature structure from the base of the wind at 10*^ g cm ^ to 10 ^ g cm ^ due to 
photoelectric heating. 
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Fig. 16. — Ionization structure of Fc from a model with LTE metal line-blanketing (top) and a model with partial non-LTE metal line- 
blanketing, including Fe I, Fe II, and Fc III (bottom). A comparison of the two structures shows they arc in good agreement in the hydrostatic 
layers where LTE is a good assumption. Further out in the wind, the two structures differ considerably. The dominant stage of Fe in the 
wind is Fe+, when treated in LTE. In non-LTE, Fe^"^ is the dominant stage of Fe throughout most of the wind, rapidly recombining to Fe*^ 
in the outermost zones. The different ionization structures in the LTE and non-LTE cases lead to significantly different synthetic line spectra 
in the UV. 
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Fig. 17. Computed number densities for specific Fe II levels from terms a^F and a'*D as a function of depth in the atmosphere. These 
levels arc the lower levels of transitions which arc found in the 1640 A to 1670 A spectral region shown in Figures 14 and 18. Level number 
densities arc shown for a model with LTE Fcl-Ill (dotted-lines) and a model with non-LTE Fel-Ill model (solid-lines). At depths greater than 
than 1 g cm~^, the hydrostatic portion of the atmosphere, the LTE and non-LTE level populuations are in close agreement. In the wind, 
non-LTE level densities are strongly depleted relative to LTE level densities by the effects of the ionizing radiation field. 
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Fig. 18. — (Top) Synthetic spectra from the two models shown in Figure 15. One model has non-LTE Fel-Ill line-blanketing, the other 
does not. The two models have otherwise identical parameters. The reduced Fe II column density in the non-LTE case means the strong 
Fe II lines probe lower velocities in the wind compared to the LTE case. For example, the two strong aborption features between 1640 A and 
1645 A are shifted by -93 kms~^ in the LTE case relative to the non-LTE case. Most strikingly, the opacity window at 1656 A is completely 
missing in the LTE case. (Bottom) Non-LTE synthetic spectrum from the top panel is compared with the lUE spectrum SWP 10375. 
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Fig. 19. — lUE spectrum LWR 07760 in a region siiowing the Mg II h and k resonance lines. Shown for comparison are synthetic spectra 
from two models: LTE Mg II and non-LTE Mg II. With the exception of the non-LTE metal line-blanketing, these models have identical 
parameters. The zero-point of the velocity scale is the h line, Aair = 2795.523 A shifted to radial velocity of Deneb, RV = —4.5 kms~^ (Evans 
1967). The interstellar components at ~ +70 kms~^, have not been modeled. 
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kms^^ which are inconsistant with the blue edges of the Mg 11 h Sz k profiles (see Figure 19). In the latter work, the 
critical assumption of a fully ionized wind is made in addition to the assumption of LTE and the Sobolev approximation. 
As shown above (§5.3, Figure 11), a fully ionized wind is inconsistent with the observed slope of the millimeter-radio 
continuum. 

We have attempted to fit the Ha profile from the HEROS red spectrum by fixing v^o at 225 kms^-'^ and computing the 
synthetic profiles in the co-moving frame from the set of radiative equilibrium expanding atmosphere models computed for 
the SED analysis. We are unable to improve upon earlier attempts to fit the observed Ha profile by adjusting any of the 
model parameters (see Figure 20). The most significant and persistent discrepancy between the synthetic and observed 
profiles is the depth of the absorption component, which is significantly weaker in the observed spectrum (residual intensity 
~ 0.6). Observed Ha profiles from 1991 (Kaufc;r et al. 1996) have a broader absorption component than shown in Figure 
20, but they are not deeper. Furthermore, while the velocity of absorption component minimum is quite variable, it is 
rarely, if ever, shifted blueward of -50 kms~^, ~20% of the terminal velocity. Models with mass-loss rates M > 10~^ 
Mq yr^^ produce synthetic profiles with absorption minima shifted to 80 kms"-'^ and beyond. Another problem is 
that in this same mass-loss range the synthetic profiles show the broad electron scattering wings prominent in other hot 
supergiants, but lacking in Deneb. Since we are unable to achieve a reasonable fit to the observed Ha profile, it is near 
impossible to reliably estimate a mass-loss rate from this line. The position of the absorption component relative to the 
models (Figure 20a) suggests a mass-loss rate in the range: 5 x 10^^ Mq yr~^ > M > 10^^ Mq yr~^. 

Higher members of the Balmcr series (H/3, H7, H5) are a better match to our synthetic profiles than Ha (see Figure 21). 
The depth of the observed line cores are well matched by our wind model synthetic profiles which have weaker cores than 
hydrostatic model synthetic profiles. However, for these same lines the wings are generally fit better by the hydrostatic 
profiles. The wind profiles also fail to match the weak metal lines in the wings of H/3 and H7. In the wind models the 
column density through the wind is clearly too large and metal lines which should form deeper in the atmosphere are 
washed out by the velocity field. The data do show the presence of a wind in the H/J and Br7 lines. The H/J line core 
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Fig. 20. — HEROS red spectrum (thick solid line) at the Ha line shown with synthetic profiles from a grid of model atmospheres with 
different parameterizations about the base model (T* = 9000 K , log3(_R*) = 1.3, M = IQ-^ Mq yr'^, (3 = 3.0, it = 15 kms"!). The 
synthetic profile for zero mass-loss is from the hydrostatic model = 8500 K , logg(i?*) = 1.3, it = 15 kms^-'. Shown are the effects on the 
synthetic profile due to the (a) luass-loss rate, (b) velocity-law exponent, (c) microturbulcnce parameter, and (d) the input radius. Positions 
of telluric lines are mark by ® symbols. All synthetic profiles have been shifted from vacuum to air wavelengths and rotationally broadened 
for vsini = 25 kms"'^. 
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is shifted sUghtly blueward from zero velocity and the Br7 profile (see Figure 21d) is noticeably asymmetric in the red 
wing, suggesting a very weak P-Cygni emission component to this line. 

For still higher members of the Balmer series, 2^17 to 2^30, the observed lines are closely matched by the wind 
model synthetic profiles (see Figure 22a). The hydrostatic model line cores are again too strong. For high members of 
the Paschen scries, 3^15 to 3^22, the difference between the wind and hydrostatic synthetic profiles is less distinct and 
both predict profiles which are a bit too strong (see Figure 22b). For the last few members of the Pfund series, 5— > 22 to 
roughly 5— >26, the wind and hydrostatic synthetic profiles are nearly indentical and are much stronger than the observed 
profiles (see Figure 22c). It is expected that for progressively weaker hydrogen series the differences between the wind 
and hydrostatic models should disappear since these weak lines form deep in the photosphere. It also appears clear that 
the n = 5 level of hydrogen in the models is overpopulated relative to Deneb's photosphere. 

The last set of line profiles we show are those from Ca II: the H ^ K lines and two members of the infrared triplet 
(IRT) (see Figure 23). The most striking characteristic of the synthetic profiles for the H, K, and A8542 lines is the 
importance of non-LTE. In LTE, Ca"*" is the dominant stage of Ca in the wind and the strong synthetic P-Cygni profiles 
of the H &: K arc evidence for this. In the non-LTE models, and apparently in Deneb's expanding atmosphere, most of 
the Ca is in the form of Ca"'""'", which has its first excited level 25 eV above the ground state, with no strong lines. For 
the H h K lines, the hydrostatic model synthetic profiles are a bit too strong and the non-LTE wind profiles are a bit 
too weak. For the model with M = 10~^ yr~^, the synthetic profiles are also significantly blueshifted, suggesting 
M< 10~^ yr~^. For the two IRT lines, only the stronger A8542 line shows a significant departure from LTE. The 
non-LTE wind model lines are a bit too strong and the hydrostatic profiles are stronger still. 

7. SUMMARY AND CONCLUSIONS 

A stellar wind module has been developed for the PHOENIX stellar atmosphere code for the purpose of computing 
atmospheric structures and detailed synthetic spectra of hot luminous stars. These models incorporate a solution of 
the expanding spherically symmetric radiative transfer equation in the co-moving frame, a treatment of detailed metal 




Fig. 21. — HEROS blue spectrum (thick solid line) at (a) H/3, (b) H7, (c) H<5, and the FTS spectrum (thick solid line with diamonds) 
at (d) Br7. Shown for comparison are synthetic profiles from three model atmospheres; two wind models, T^, = 9000 K, log (/(/?*) = 1.3, 
f) = 3.0, it = 15 kms~l, M = 10~^ M© yr^^ (dotted line) and M = 10^"^ M© yr-l(dashed line), and the hydrostatic model (dash-dot 
line) T^, = 8500 K , \ogg{R^,) = 1.3, it = 15 kms'^ The synthetic profiles for comparison with the HEROS data have been shifted from 
vacuum to air wavelengths. All synthetic profiles have been rotationally broadened for vsini = 25 kms~^. 
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Fig. 22. — Observed hydrogen lines (thick soUd lines), from three scries: (a) Balrncr, (b) Paschen, and (c) Pfund. These lines are all high 
members of there respective series: (a) 2^17 to 2^30, (b) 3^15 to 3^22, and (c) 5^22 to ~-5^26. Other strong features in the spectra 
include Or II A3678, Ti II A3685, O I A8446, and Ca II AA8498,8542. The Paschen series is contaniinatcd by many weak telluric features. The 
observed Balmer and Paschen series are from the HEROS blue and red spectra respectively and the Pfund series is from the FTS spectrum. 
Shown for comparison are synthetic profiles from two model atmospheres: a wind model (dashed line), = 9000 K, logg(i?*) = 1.3, /3 = 3.0, 
= 15 kms"\ M = 10"*' M0 yr^^ and a hydrostatic model (dotted hne) = 8500 K , logg(R*) = 1.3, = 15 kms"^. The synthetic 
profiles for comparison with the HEROS data have been shifted from vacuum to air wavelengths. All synthetic profiles have been rotationally 
broadened for vsini = 25 kms~^. 
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Fig. 23. — HEROS spectra (thick solid lines) at the Ca II lines (a) K, (b) H, (c) A8498, and (d) A8542. Shown for comparison are synthetic 
profiles from four model atmospheres: two wind models with non-LTE Cal-Ill, plus one wind model and one hydrostatic model with LTE 
Cal-Ill. The wind models: T^, = 9000 K, logg(i?*) = 1.3, /3 = 3.0, = 15 kms^^, M = 10"^' Mq yr"! (long dash-dot fine) and M = lO"'' 
Mq yr"^(LTE dashed line, non-LTE dotted line) and the hydrostatic model: T* = 8500 K , logg(ii*) = 1.3, $t = 15 kms"^ (short dash-dot 
line). The synthetic profiles have been shifted from vacuum to air wavelengths. All synthetic profiles have been rotationally broadened for 
Dsinj = 25 kms~^. 
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line-blankcting in non-local thermodynamic equilibrium, and a convergence of the model temperature structure under the 
constraint of radiative equilibrium. We have used the model atmospheres and synthetic spectra from this code to better 
illuminate what we do not understand about the spectral energy distribution, the mass-loss rate, and the fundamental 
stellar parameters of the A-type supcrgiant Deneb. 

7.1. Millimeter and Centim,eter Radiation Deteetions 

We have reported the first detections of Deneb at millimeter and centimeter wavelengths with the SCUBA and the 
VLA. These new data, in conjunction with IR photometery from 2 /xm to 60 /xm, provide a better glimpse at the thermal 
continuum from Deneb's extended atmosphere. The slope of the radio spectrum implies a stellar wind where hydrogen is 
partially ionized, a result reproduced by our models. We find that the effects of line-blanketing on the synthetic millimeter 
and radio continua are significant and must be taken into when analyzing the continua of partially ionized winds. 

7.2. Angular Diameter Measurement and Limb-darkening 

Visibility measurements from the NPOI indicate a uniform-disk angular diameter for Deneb of ^ud = 2.40 ± 0.06 mas. 
Our models predict that the center-to-limb profile for Deneb in the /-band is significantly extended relative to hydrostatic 
model predictions. We believe this limb profile causes the uniform disk angular diameter to be larger than the angular 
diameter ^ross = 2.32 ± 0.06 mas corresponding to i?Ross) the radius at the Rosseland mean optical depth trqss = 2/3. 
The result is a correction to the uniform disk diameter which is less than unity, contrary to the limb- darkening correction 
made in the calculation of the effective temperature for more compact stars. This limb profile prediction will hopefully 
be testable in the coming years as longer base-line optical/IR interferometers come on-line. 

7.3. Interferorn,etie Effective Temperature 

Considering the uncertainites in the spectrophotometry, the interstellar extinction, and the corrected NPOI angular 
diameter, we find an inteferometric fundamental effective temperature of T^^^ = 8600 ± 500 K for Deneb. The biggest 
uncertainty in the value of the effective temperature is the uncertainty in the color excess, E{B — V^) = O.OOi'^Q Qg. This 
effective temperature together with the distance estimates from Hipparcos and Deneb's Cyg 0B7 association membership, 
and theoretical evolutionary tracks yield: i?Ross — 180 Rq, i ~ 1.6 ± 0.4 x 10* Lq, M ~ 20 — 25 Mq, and log^ ~ 1.3. 

7.4. Least-squares Spectrum Analysis 

A least-squares comparison of 55 synthetic spectral energy distributions from 1220 A to 3.6 cm with the spectral energy 

distibution of Deneb provide estimates (Icr) for the star's effective temperature and mass- loss rate of: T^i°^^ — 8420 ±100 

K and M ~ 8 ± 3 x 10"'' Mq yr~^. The best fitting models to the SED indicate that the color excess toward Deneb is 
E{B — V) = 0.06±0.01. This effective temperature is consistent with the inteferometric fundamental effective temperature. 
The smaller la error in the least-squares value of T^^^ relative to interferometric value reflects the tighter contraint on 
the color excess from the least-squares fit. However, both values of T^^^ are subject to systematic errors due to the 
poorly established interstellar extinction curve toward Deneb. 

The least-squares analysis shows that, in addition to the effective temperature and mass-loss rate, the steepness of 
/3-law velocity field, the microturbulence, and the surface gravity all have a significant effect on the model temperature 
structure and ultraviolet spectral energy distribution, ft is interesting to point out that models with /? outside the range 
3.0 ± 1.0 provide very poor fits to the UV spectrophotometry and that /3 values near 3.0 are not consistent with the 
radiation driven wind models of Achmad et al. (1997), which predict /? < 1 for A-type supergiant winds. 

7.5. Non-LTE Effects on the Line Spectrum 

The mass-loss rate range derived from the least-squares analysis is reasonably consistent with that derived from high 

dispersion spectra of Deneb when non-LTE metal- line blanketing is treated, in particular for species Fei-iii, Mgii, and 
Cai-iii. The primary effect of the non-LTE treatment is to enhance the ionization of these elements in the wind relative 
to LTE. The non-LTE model ionization structure produces much lower column densities of Fe II, Mg II, and Ca II relative 
to LTE and causes strong lines to form deeper in the wind at lower velocities, in better agreement with the observations. 
For the Fe II UV lines and the Ca II H K lines, the non-LTE effects are striking. The mass-loss rate inferred from the 
Fe II UV lines is 50 to 100 times larger in non-LTE than in LTE. The synthetic H Sz K profiles show strong synthetic 
P-Cygni profiles in LTE which are not present in the data or in the non-LTE synthetic profiles. 

7.6. Problems with Ha and Other Hydrogen Lines 

We are unable achieve a reasonable fit to a typical Ha P-Cygni profile for Deneb with any model parameters over a 
reasonable range. The most prominent problem is that the model absorption component is always too strong. While 
Deneb's Ha profile is variable, it is not hugely variable, and the model Ha profiles don't look like the observed profile at 
anytime. The wind model line profiles generally provide a better match to the higher Balmer series lines than hydrostatic 
model line profiles. Both wind models and hydrostatic models match the higher Paschen series lines equally well. Both 
types of models fail badly for the higher Pfund series lines. These failures indicate that a spherically symmetric, expanding, 
steady state, line-blanketed, radiative equilibrium structure is not consistent with the conditions under which a typical 
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Ha and the higher Pfund Unes form. The strong hncs of Fe II, Mg II, and Ca II that we can reproduce with our models 
originate from the ground state or low-exitation levels and apparently are not sensitive to the conditions which affect the 
2^3 transition of hydrogen. Our inability to fit a typical Ha profile from Deneb is very bothersome because Ha is the 
most commonly used mass-loss rate diagnositic for A-type supergiants. 

7.7. Mass-loss Rate 

The thermal millimeter and radio continua of Deneb when compared with our models provide both lower and upper 
limits to its mass-loss rate. The detection of radiation at 3.6 cm, far in excess of what is expected for no mass-loss, 
provides the lower limit at about 10"'' Mq yr~^. The 3a upper limit to the 870 fxm flux provides the upper limit at about 
10~^ Mq yr~^. If Deneb's extended atmosphere departs significantly from spherical symmetry, this would systematically 

bias our estimates for the mass-loss rate. However, there is no evidence for deviations from spherical symmetry from 
instrinsic polarization studies, either in the continuum or in the spectral lines H/3, H7, Ca II H (Clarke & Brooks 1984, 
and references therein). 

The line spectrum of Deneb provides mainly upper limits to the mass-loss rate. While the observed P-Cygni character 
of the Ha and Mg II h and k lines are clear indicators of the stellar outflow, a hydrostatic model with no mass loss is 
really not that bad a fit to most of the UV line spectrum. It is the more subtle characteristics of the UV lines, their 
widths and desaturated profiles which show evidence of the velocity field. Many weak metal lines in the optical spectrum 
are also very well matched by models with no mass-loss and all our models with M> 10~^ M0 yr~^ predict profiles for 
these lines which on average tend to be washed out and too shallow. Certainly the column densities in many of these 
transitions are too high. While non-LTE ionization effects solve this problem for some lines, others are still problematic. 
A line-by-line study with a complete grid of models which include Til-iii, Cri-iii, and other metals in non-LTE is needed 
to analyze this issue fully. We have model atoms for these species and will investigate this in future work. 

While previously published estimates for Deneb's mass-loss rate range over more than three orders of magnitude, we 
can with confidence reduce this uncertainty to one order of magnitude, but not much better. Reducing this uncertainty 
further requires a better match simultaneously to the spectral energy distribution and the spectral lines. This may require 
removing several of our simplifying assumptions such as radial symmetry, homogeniety (filling- factor equal to unity), 
time-independence, and no mechanical dissipation. 

We suggest that before the Wind Momentum-Luminosity Relationship, thus far based solely on Balmer line fits, can be 
firmly established for A-type supergiants, researchers should work to toward firmly establishing fundamental parameters 
and stellar wind properties of Galactic A-type supergiants by checking that models are consistant with both the spectral 
energy distributions and the line spectra. 
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APPENDIX 

A. CONSTRUCTION OF MODEL ATMOSPHERES 

The stellar wind module developed by Aufdenberg (2000) for the PHOENIX code links a unified model atmosphere 
structure to the non-LTE line-blanketing and expanding atmosphere radiative transfer modules already in the code. 
Important improvements over previous models lie in the solution of the radiative transfer equation, the degree of consistent 
LTE and non-LTE metal line-blanketing, and the solution of the temperature structure from the condition of energy 
conservation. 
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A.l. Line-blanketing Specifics and Non-LTE Model Atoms 

All models presented in this paper include in LTE approximately 7 x 10^ of the most important blanketing lines 
dynamically selected from a list of 42 x 10*' lines (Kurucz & Bell 1995) from 39 elements with up to 26 ionization stages 
in addition to the strongest 517 lines for the spectra H I, He I, He II, which are treated in non-LTE. More sophisticated 
models include an additional 2.3 x 10* lines in non-LTE for these spectra: C I, C II, C III, N I, N II, N III, O I, O 11,0 III, 
Mg I,Mg II, Mg III, Ca I, Ca II, Ca III, Fc I, Fc II, Fc III. The number of levels and lines from each species are shown in 
Table A8. For the complete set of non-LTE model ions in the PHOENIX code see Short ct al. (1999). The references to the 
atomic data used to construct the model atoms and ions are given in Table A9. 

The non-LTE lines from the model atoms and ions replace the LTE lines and all the lines that are not in non-LTE are 
handled in LTE as an approximation. Details on non-LTE blanketing techniques can be found in (Short et al. 1999, and 
references therein). Even in our most sophisticated models the vast majority of lines arc treated in LTE. The bulk of the 
total opacity is, however, computed in detailed non-LTE. The line-blanketing is true linc-blankcting in the sense that the 
effect of these millions of lines on the temperature structure is computed through the condition of radiative equilibrium 
(see below). Therefore, the line-blanketing is consistently included in both the model atmosphere and final synthetic 
spectrum calculations. 

A.2. The Wind Module 

We treat the problem of radiation-driven winds with the following simplifying assumptions: (1) time- independent, 
steady-state; (2) no magnetic fields or rotation; (3) radial symmetry; and (4) a smooth single-phase flow with no mechanical 
dissipation, heat conduction, viscosity, or advection. 

A.2.1. Atmospheric Structure and Boundary Conditions 

The gas escaping from the model star is treated as a hydrodynamic flow, satisfying the conservation of mass, momentum, 
and energy. The mass conservation is specified by the equation of continuity, 

M = Anr'^ p{r)v{r) (Al) 

where M is the steady state mass-loss rate, p{r) is the gas density, and v{r) is the velocity field. The velocity field is 
specified by a /3-law of the form 

v{r)=Voo{l-R./rf. (A2) 

The layers of the atmosphere below i?* are modeled as hydrostatic and in these layers the velocity is set to zero. 

The structure of the outer wind layers is computed on a fixed radial grid. The structure of the inner hydrostatic layers 
of the atmosphere is computed on an optical depth grid. Typically 50 layers are used consisting of 35 layers for the wind 
region and 15 layers for the hydrostatic region. The radial grid points are prescribed by a power-law in order to finely 
sample the inner region of the wind (Steffen et al. 1997), where the velocity gradient is steepest, and to coarsely sample 
the outer portion of the wind where the velocity gradient is small. 

Table A8 

Number of Le\"els and Primaky Transitions''' Treated in Non-LTE 



Stage of Ionization 



Element 


I 


II 


III 


H 


30/435 






He 


19/37 


10/45 




C 


228/1387 


85/336 


79/365 


N 


252/2313 


152/1110 


87/226 


O 


36/66 


171/1304 


137/765 


Mg 


273/835 


72/340 


91/656 


Ca 


194/1029 


87/455 


150/1661 




li)i/()i)o:! 


()i7/i:!()75 


.")()() /I)72i 



■^Primary transitions are those that connect ob- 
served levels and have log(g'/) values greater than 
-3. 

''Values for the complete models of stages i-iii 
of Fe are shown. In our calculations, only the first 
200 levels of each stage are used. 
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A. 2. 2. Wind region 
The radii r(Z) in layers I = 1, are specified by, 

r(Z)=r(l) + {r(U-r(l)} 



vC-i) 



(U-i) _ 1 



(A3) 



where is the layer (typically layer 35) immediately above the dynamic-static transition. Layer 1 is the outermost 

layer and typically wc choose the radial extension in the wind to be a factor of 200, r(l) = 200i?*. The value of r(Z^,) 
is itcratively selected to provide a predetermined degree of smoothness in the density profile across the dynamic-static 
transition. Generally the jump in density must be less than a factor of 5 across the transition. The input parameter a 
may be adjusted to optimize the distribution of radial points for a specific velocity field (typically, a ~ 1.5). 

The density p{r) in the expanding layers is specified from the mass flux continuity equation (Al), where the velocity 
field has the form of a /3-law from equation (A2). The input radius ii* and the model effective temperature Tj, set the 
luminosity L of the star, 

L = A-nRlaT^. (A4) 

Since we use a form of the /3-law lacking a macroturbulent velocity, the switch in the structure setup from the dynamic 

layers to the static layers occurs near the point in the atmosphere where the wind velocity, v{r), drops below the micro- 
turbulent velocity, An initial guess at this radius, r(/*) > i?*, is calculated by finding the radius at which v{r) = 
from equation (A2): 

K^*) = 7 (A5) 



This radius is then adjusted inward slightly to improve the degree of smoothness of the density profile across the dynamic- 
static transition. In practice, this results in a discontinuity in the velocity field across dynamic-static boundary which is 
very small (~0.2 kms~^), much smaller than the microturbulent width of the lines, > 5 kms~^. 



In the hydrostatic layers {I 
grid. 



A. 2. 3. Hydro st,at,ic region 
li, + 1, ...,N), the structure is computed on a fixed logarithmically spaced optical depth 



log{r(0} = log{r(/.)} + (Z - Q 



/ log{T(iV)}-log{T(Z.)} 



(A6) 



The inner continuum extinction optical depth boundary, t{N) = r(2000A) = 10^, is set well below the thermalization 
depth. The pressure structure is computed by numerically integrating the hydrostatic equation, 

dP _ ffeff 

dr K 



(A7) 



Table A9 

Atomic Data Sources for Model Atoms and Ions'' 



Model Species 


Levels 


b — h cross-sections 


b — f cross-sections 






Radiative Collisional^ 


Photoionization^ Collisional 



H I 


Johnson 


Johnson 


Johnson 


M84 


Drawin 


He I 


Martin 


WSG 


BK, Lanzafame 


M84 


MS 


He II 


Johnson 


Johnson 


Johnson 


M84 


MS 


C i-iii 


KB 


KB 


Allen, VR 


BRP, RM 


Drawin 


N i-iii 


KB 


KB 


Allen, VR 


BRP, RM 


Darwin 


O I 


KB 


KB 


Allen, VR 


BRP, RM 


CLL 


O ii-iii 


KB 


KB 


Allen, VR 


BRP, RM 


Drawin 


Mg I-III 


KB 


KB 


Allen, VR 


RM 


Drawin 


Ca I-III 


KB 


KB 


Allen, VR 


BRP, RM 


Drawin 


Fe I-III 


KB 


KB 


Allen, VR 


BRP, RM 


Drawin 



'References - Alien = Allen (1973); BK = Bcrrington & Kingston (1987); BRP = Bautista ct al. 
(1998); KB = Kurucz & Bell (1995); CLL=Chung et al. (1993, 1994); Drawin = Drawin (1961); 
Lanzafame = Lanzafame et al. (1993); M84 = Mathisen (1984); MS = Mihalas & Stone (1968); 
Johnson = Johnson (1972); RM = Reilman &: Manson (1979); VR = Van Rcgcmortor (1962); WSG 
= Wioso ot al. (1966) 

^Reference to .Allen is only for cross-sections between bound levels not coupled by radiative 
transitions 

^Reference to BRP is only for ground-state cross-sections 
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where ^eff is the effective gravitational acceleration 

fifeflF = 9{R*)^ - 9T^d{r) (A8) 

the difference of the Newtonian acceleration and the radiative acceleration grad computed from 

1 

where XA(»")ffA(T') is the product of the total line and continuum opacity with the first moment of the radiation field 
provided by solution of the radiative transfer equation (sec below). 
The radii in the hydrostatic layers are computed by integrating 

dr = — dr. (AlO) 

Once the velocity field and radii are computed for all layers, the non-grey, non-LTE radiative transfer problem is solved 
and the temperature structure is adjusted iterativcly to achieve energy conservation. 

A. 3. Radiative Transfer and Radiative Equilibrium 

The spherically symmetric radiative transfer equation (SSRTE) in the co-moving (Lagrangian) frame (Mihalas & Weibel- 
Mihalas 1984) 



A/ 



+7 {^^^±4^ + 7^(1 + + 2/3M)f } / 



(All) 



is solved using an Operator Splitting/ Approximate A-opcrator Iteration (OS/ALI) method using short-characteristics 
(Hauschildt 1992a). Improvements since this work include a more accurate discretization of the d{XI)/d\ term needed to 
improve the numerical accuracy necessary for /3-law winds. 

In equation (All), /i is the direction-cosine, / is the specific intensity, the velocity is measured in units of the speed of 
light c, (3{r) = v{r)/c, and 7 = 1/(1 — /3^)^/^. The principal reason for choosing to solve the radiative transfer problem 
in the co-moving, fluid rest frame is that although the frequency spectrum may be complex, with many overlapping 
lines, the opacity is isotropic. This leads to great simplification in the radiation-matter interaction terms relative to the 
observer's frame equation. For example, the flow velocity is not an argument of the line profile function in the co-moving 
frame. Frequencies are measured in the frame of the moving fluid, therefore wc can use a static line profile. In this frame 
one can directly integrate over angle and use co-moving-frame moments of the radiation field in the computation of the 
radiation pressure and radiative equilibrium equations. The angular dependence of the opacity in the observer's frame 
leads to computational difficulties, and in flows with small velocity gradients and many overlapping spectral lines, escape 
probability methods, like the Sobolev approximation, are not valid. 

Wc iterativcly compute the temperature structure from the condition of energy conservation in the co-moving frame. 
The conservation of energy, assuming radiative equilibrium, can be specified in the Lagrangian frame (Hauschildt 1992b) 

by 

poo 

{vx-XxJx)dX + S = (A12) 





where rjx and Jx describe the emissivity and mean intensity respectively. We include the kinetic and gravitational energy 
of the wind, S, in the energy equation, although they are less than 0.1% of the total nuclear energy generated by the star 
for the models presented here and thus have negligible effect on the temperature structure. For very high wind momenta 
and small luminosities, the energy required to drive the wind might be large enough to affect the temperature structure. 
Our numerical solution of the special relativistic radiative energy equation is described in Hauschildt & Baron (1999). 

Energy conservation errors are typically < 1% in both the luminosity and luminosity derivative in each layer. Models 
are started from a grey temperature structure and iterations are performed with continuum opacity only. Next, more 
iterations are performed with the addition of full LTE metal-line blanketing opacity. This yields a converged LTE linc- 
blanketed model. Non-LTE models are converged starting from converged LTE line-blanketed models. Non-LTE species 
are added and the model structure is converged in a series of steps: first H and He; next C, N, O; next Mg, Ca, Fe, etc. 
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